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Forward production of charged pions with incident protons on nuclear 

targets at the CERN PS 



Measurements of the double-differential tt^ production cross-section in the range of 
momentum 0.5 GeV/c < p < 8.0 GeV/c and angle 0.025 rad < 9 < 0.25 rad in collisions 
of protons on beryllium, carbon, nitrogen, oxygen, aluminium, copper, tin, tantalum and 
lead are presented. The data were taken with the large acceptance HARP detector in 
the T9 beam line of the CERN PS. Incident particles were identified by an elaborate 
system of beam detectors. Thin targets of 5% of a nuclear interaction length were used. 
The tracking and identification of the produced particles were performed using the 
forward system of the HARP experiment. Results are obtained for the 
double-differential cross-sections d 2 a/dpdQ mainly at four incident proton beam 
momenta (3 GeV/c, 5 GeV/c, 8 GeV/c and 12 GeV/c). Measurements are compared 
with the GEANT4 and MARS Monte Carlo generators. A global parametrization is 
provided as an approximation of all the collected datasets which can serve as a tool for 

quick yields estimates. 
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PACS numbers: 



I. INTRODUCTION 



Final HARP measurements of the double- 
differential cross-section, d 2 er/dpdSl, for 7r ± for- 
ward production (in the range of momentum 
0.5 GeV/c < p < 8.0 GeV/c and angle 0.025 rad < 
9 < 0.25 rad) by incident protons of 3 GeV/c, 
5 GeV/c, 8 GeV/c, 8.9 GeV/c (Be only), 12 GeV/c 
and 12.9 GeV/c (Al only) momentum impinging on 
thin solid beryllium, carbon, aluminium, copper, 
tin, tantalum or lead targets of 5% nuclear inter- 
action length are presented. Our results on the for- 
ward production of ir + in p-Al interactions at 12.9 
GeV/c and p-Be interactions at 8.9 GeV/c 0, 
useful for the understanding of the the K2K (3|] and 
MiniBoone, SciBooNE neutrino beams, have been 
reanalyzed now with a consistent binning. In ad- 
dition, results at 12 GeV/c on thin carbon targets 
and on cryogenic N2 and O2 targets, relevant for a 
precise calculation of atmospheric neutrino fluxes 
and improved modeling of extensive air showers, 
previously reported in references [H and are 
included for completeness. 

The HARP experiment at the CERN PS was 
designed to measure hadron yields from a large 
range of solid and cryogenic nuclear targets and 
for incident particle momenta from 1.5 GeV/c to 
15 GeV/c. This corresponds to a proton momen- 
tum region of great interest for neutrino beams and 
far from coverage by earlier dedicated hadropro- 
duction experiments @, Q ■ The main motivations 
were to measure pion yields to allow substantially 
improved calculations of the atmospheric neutrino 
flux Q to be made, to measure particle yields as in- 
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put for the flux calculation of accelerator neutrino 
experiments [10( , to help to design the targetry for 
a future neutrino factory [ll[ and to improve the 
reliability of extensive air shower simulations, by 
reducing the uncertainties of hadronic interaction 
models in the low energy range [12j . 

Covering an extended range of solid targets in the 
same experiment, it is possible to perform sys- 
tematic comparisons of hadron production mod- 
els with measurements at different incoming beam 
momenta over a large range of target atomic num- 
ber A. Pion production data at low momenta 
(< 25 GeV/c) are extremely scarce and HARP 
is the first experiment to provide a large data 
set, taken with many different targets, full par- 
ticle identification and large detector acceptance. 
These data, together with the ones published in 
reference (l3j . are also of great interest for correc- 
tions of secondary interactions in detector studies 
in particle physics experiments. 

Existing data are mainly at fixed production an- 
gles with Be targets [3] and are affected usuall y by 
large uncertainties. Only two experiments (TEL [la] 
at 19.2 GeV/c and 24 GeV/c were performed with 
an extended set of nuclear targets. The E910 ex- 
periment at BNL has recently published data at 
12.3 GeV/c and 17.5 GeV/c with protons on Be, 
Cu and Au targets [T3|. 

Data were taken in the T9 beam of the CERN PS. 



II. EXPERIMENTAL APPARATUS 



The HARP experiment makes use of a large- 
acceptance spectrometer consisting of a forward 
and a large-angle detection system. The HARP 
detector is shown in Fig. [T] and is fully described 
in Ref. [T8l ]. The forward spectrometer - based on 
five modules of large area drift chambers (NDC1- 
5) [Til and a dipole magnet complemented by a 
set of detectors for particle identification (PID): a 
time-of-fiight wall (TOFW) [H, a large Cherenkov 
detector (CHE) and an electromagnetic calorime- 
ter (ECAL) - covers polar angles up to 250 mrad, 
which is well matched to the angular range of in- 
terest for the measurement of hadron production 
to calculate the properties of conventional neutrino 
beams. As the most downstream detection element 
visible in the figure, the 1.4 m wide beam muon 
identifier (BMI) is used to measure the muon con- 
tamination in the beam. 
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The discrimination power of TOFW time-of-flight 
below 3 GeV/c and the Cherenkov dectector above 
3 GeV/c are combined to provide powerful separa- 
tion of forward pions and protons. The calorime- 
ter is used only for separating pions and electrons 
when characterizing the response of the other de- 
tectors. The resulting pion identification efficiency 
is bigger than 96% (98%) for pion momenta larger 
than 0.5 (3.0) GeV/c, with purity around 99.5%. 

The large-angle spectrometer - based on a Time 
Projection Chamber (TPC) and Resistive Plate 
Chambers (RPCs) located inside a solenoidal mag- 
net - has a large acceptance in the momentum and 
angular range for the pions relevant to the produc- 
tion of the muons in a neutrino factory. For the 
analysis described here only the forward spectrom- 
eter and the beam instrumentation are used. 

The HARP experiment took data in 2001 and 
2002. The momentum definition of the T9 beam 
is known with a precision of the order of 1% [2lj ]. 

The target is placed inside the inner field cage of 
the TPC, in an assembly that can be moved in 
and out of the solenoid magnet. The solid targets 
used for the measurements reported here have a 
cylindrical shape with a nominal diameter of about 
30 mm. Their thickness along the beam is equiva- 
lent to about 5% of an interaction length Ai |41j . 

A set of four multi-wire proportional chambers 
(MWPCs) measures the position and direction of 
the incoming beam particles with an accuracy of 
«1 mm in position and «0.2 mrad in angle per 
projection. A beam time-of-flight system (BTOF) 
measures the time difference of particles over a 
21.4 m path-length. It is made of two identi- 
cal scintillation hodoscopes, TOFA and TOFB 
(originally built for the NA52 experiment [22]), 
which, together with a small target-defining trigger 
counter, also used for the trigger, provide particle 
identification at low energies. This provides sepa- 
ration of pions, kaons and protons up to 5 GeV/c 
and determines the initial time at the interaction 
vertex (to). The timing resolution of the combined 
BTOF system is about 70 ps. A system of two 
N 2 -filled Cherenkov detectors (BCA and BCB) is 
used to tag electrons at low energies and pions at 
higher energies. The electron and pion tagging effi- 
ciency is found to be close to 100%. A set of trigger 
detectors completes the beam instrumentation. 

The positive beam used in this analysis contains 
mainly positrons, pions and protons, with a small 
admixture of kaons and deuterons and heavier ions. 
The proton fraction in the incoming beam goes 
from - 35% at 3 GeV/c to - 92% at 12 GeV/c. 
The selection of beam protons is done requiring a 
pulse height consistent with the absence of signal in 



both beam Cherenkov detectors (BCA and BCB). 
At the lowest beam energy, 3 GeV/c, the BTOF is 
used to reject pions which at that momentum do 
not produce Cherenkov light. At 5 GeV/c, the ir/p 
separation is obtained by the BTOF system and 
only one of the Cherenkov (usually BCB), while 
the other is used to tag e ± . Ions are rejected by the 
BTOF system at all momenta. In the 12.9 GeV/c 
beam the two Cherenkov counters were operated 
with different pressures to make it possible to mea- 
sure the kaons separately from the protons and pi- 
ons. After tagging the beam contamination from 
all particle species is negligible. Only events with 
a single reconstructed beam track in the four MW- 
PCs, good timing measurements in BTOF and no 
signal in the beam halo counters are accepted. 

A downstream trigger in the forward scintillator 
trigger plane was required to record the event, ac- 
cepting only tracks with a trajectory outside the 
central hole (60 mm diameter) which allows beam 
particles to pass. Accepting only tracks with a 
trajectory outside the central hole, the efficiency 
is measured to be >99.8%. Particle identification 
is done later, at the analysis stage, via the down- 
stream PID detectors. 

The length of the accelerator spill is 400 ms with 
a typical intensity of 15 000 beam particles per 
spill, and less for the lower momentum settings. 
The average number of events recorded by the data 
acquisition ranges from 300 to 350 per spill. 

The absolute normalization of the cross-section 
was performed using 'incident-proton' triggers. 
These are triggers where the same selection on the 
beam particle was applied but no selection on the 
interaction was performed. The rate of this trig- 
ger was down-scaled by a factor 64. These un- 
biased events are used to determine A pot in the 
cross-section formula (1), see later. 



III. DATA ANALYSIS 



The analysis procedure is similar to the one re- 
ported in Refs. 0], [|, and [I3| where with 
respect to our first paper on forward pion produc- 
tion in p-Al interactions [l| , a number of improve- 
ments to the analysis techniques and detector sim- 
ulation had been made. Tracks are reconstructed 
in the drift chambers downstream of the magnet. 
Only tracks with at least one hit in the TOFW 
are accepted for the analysis. Hits are searched for 
in the Cherenkov detector consistent with these 
tracks to complete the particle identification. Sec- 
ondary track selection criteria are optimized to en- 
sure the quality of momentum reconstruction and 
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TOFW 



TPC and RPCs in 
solenoidal magnetj 
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FIG. 1: Schematic layout of the HARP detector. The convention for the coordinate system is shown in the 
lower-right corner. The three most downstream (unlabelled) drift chamber modules are only partly equipped 
with electronics and are not used for tracking. The detector covers a total length of 13.5 m along the beam axis 
and has a maximum width of 6.5 m perpendicular to the beam. 



a clean time-of-flight measurement while maintain- 
ing a high reconstruction efficiency. The method 
to combine the information from these PID detec- 
tors is described in (23J. In the kinematic range 
of the current analysis the pion identification ef- 
ficiency is about 98%, while the background from 
mis-identified protons is well below 1%. The mo- 
mentum reconstruction is performed by means of 
a special implementation 24j of the Kalman fil- 
ter [25[ which uses the position of the target as a 
constraint. For the final stage of the analysis the 
same unfolding technique as in Refs. [f| and 
[l3| has been applied. 

The background induced by interactions of beam 
particles in the materials outside the target is mea- 
sured by taking data without a target in the target 
holder ("empty target data"). These data are sub- 
ject to the same event and track selection criteria 
as the standard data sets. To take into account 
this background the number of particles of the ob- 
served particle type in the "empty target data" 
are subtracted bin-by-bin (momentum and angu- 
lar bins) from the number of particles of the same 
type. The uncertainty induced by this method is 
labeled "empty target subtraction" . The collected 
event statistics on the different solid targets is sum- 
marised in Table [T] 



Cross-section calculation 



The double-differential cross-section is calculated 
as follows: 



d 2 a° 



(Pi, 



A 



1 



1 



dpdtt^ u 01 N A pt iV po t ApiAtij 



(1) 



where 



p'. ,0'- ,a' 



j(Pi,0j) is the cross-section 



dpdQ W« ' "3 1 

mb/(GeV/c sr) for the produced parti- 
cle type a (p, 7r + or 7r~) for each true 
momentum and angle bin (pi,6j) covered in 
this analysis; 

N a (p'i,0'j) is the number of reconstructed 
particles of type a' in bins of reconstructed 
momentum p\ and angle 6j in the raw data, 
after subtraction of empty target data (due 
to beam protons interacting in material other 
than the nuclear target). These particles 
must satisfy the event, track and PID selec- 
tion criteria. 



• ■M.p^eap'e'a' i s t ne correction matrix which ac- 
counts for the finite efficiency and resolution 
of the detector. It unfolds the true vari- 
ables pi ,0j,a from the reconstructed vari- 
ables p[, Oj,a' and corrects the observed 
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TABLE I: Total number of events and tracks used in the various nuclear 5% Ai solid target data sets and the 
number of accepted protons on target as calculated from the pre-scaled incident beam triggers. Numbers are for 
incident proton in units of 10 events. 



Data set (GeV/c) 


3 


5 


8 


8.9 


12 


12.9 


Total DAQ events 


(Be) 


1113 


1296 


1935 


5868 


1207 






(C) 


1345 


2628 


1846 




1062 






(Al) 


1159 


1789 


1707 




619 


4713 




(Cu) 


624 


2079 


2089 




745 






(Sn) 


1637 


2828 


2404 




1803 






(Ta) 


1783 


2084 


1965 




866 






(Pb) 


1911 


2111 


2266 




487 




Accepted beam protons 


(Be) 


99 


289 


761 


2103 


580 




with forward interaction (C) 


101 


542 


709 




470 






(Al) 


86 


376 


637 




306 


2116 




(Cu) 


73 


408 


741 




363 






(Sn) 


217 


528 


818 




856 






(Ta) 


281 


398 


668 




403 






(Pb) 


310 


387 


758 




221 




Final state tt (tv + ) 


(Be) 


0.08 (0.4) 


1.2 (2.9) 


8.0 (15.5) 


26.5 (48.6) 


11.0 (18.0) 




selected with PID 


(C) 


0.06 (0.3) 


1.9 (4.9) 


6.9 (13.8) 




8.2 (13.5) 






(Al) 


0.05 (0.2) 


1.3 (3.3) 


6.4 (12.9) 




5.8 (9.3) 


45.9 (70.9) 




(Cu) 


0.03 (0.1) 


1.3 (3.1) 


7.1 (12.6) 




6.8 (10.4) 






(Sn) 


0.07 (0.2) 


1.5 (3.4) 


7.6 (13.6) 




15.2 (23.5) 






(Ta) 


0.08 (0.2) 


1.0 (2.2) 


5.5 (9.7) 




7.1 (10.4) 






(Pb) 


0.08 (0.2) 


0.9 (2.0) 


5.7 (10.3) 




3.7 (5.4) 
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particle number to take into account ef- 
fects such as reconstruction efficiency, accep- 
tance, absorption, pion decay, tertiary pro- 
duction, PID efficiency and PID misidcntifi- 
cation rate. 

• A> l\rb and Ap.An, are normalization 
factors, namely: 

is the number of target nuclei per 
unit area (4^ |: 

iVp t is the number of protons on target; 

Api and Afl,- are the bin sizes in mo- 
mentum and solid angle, respec- 
tively [H[. 

We do not make a correction for the attenuation 
of the proton beam in the target, so that strictly 
speaking the cross-sections are valid for Ai = 5% 
targets. 

The calculation of the correction matrix 
Mp°$. ap >o' a ' 1S done with the unfolding method 
introduced by D'Agostini [26| Q. This method 
has been used in the recent HARP publica- 
tions [J, [H, [lH and it is also applied in the analysis 
described here. 

The Monte Carlo simulation of the HARP setup 
is based on GEANT4 [13]. The detector materials 
are accurately described in this simulation as well 
as the relevant features of the detector response 
and the digitization process. All relevant physics 
processes are considered, including multiple scat- 
tering, energy loss, absorption and re-interactions. 
The simulation is independent of the beam particle 
type because it only generates for each event ex- 
actly one secondary particle of a specific particle 
type inside the target material and propagates it 
through the complete detector. A small difference 
(at the few percent level) is observed between the 
efficiency calculated for events simulated with the 
single-particle Monte Carlo and with a simulation 
using a multi-particle hadron-production model. 
A similar difference is seen between the single- 
particle Monte Carlo and the efficiencies measured 
directly from the data. A momentum-dependent 
correction factor determined using the efficiency 
measured with the data is applied to take this 
into account. The track reconstruction used in 
this analysis and the simulation are identical to 
the ones used for the ir + production in p-Be colli- 
sions [H . A detailed description of the corrections 
and their magnitude can be found there. 

The reconstruction efficiency (inside the geometri- 
cal acceptance) is larger than 95% above 1 .5 GeV/c 
and drops to 80% at 0.5 GeV/c. The require- 
ment of a match with a TOFW hit has an effi- 
ciency between 90% and 95% independent of mo- 



mentum. The electron veto rejects about f % of the 
pions and protons below 3 GeV/c with a remaining 
background of less than 0.5%. Below Cherenkov 
threshold ( - 3 GeV/c) the TOFW separates pi- 
ons and protons with negligible background and 
an efficiency of «98% for pions. Above Cherenkov 
threshold the efficiency for pions is greater than 
99% with only 1.5% of the protons mis- identified 
as a pion. The kaon background in the pion spec- 
tra is smaller than 1%. 

The absorption and decay of particles is simulated 
by the Monte Carlo. The generated single particle 
can re-interact and produce background particles 
by hadronic or electromagnetic processes, thus giv- 
ing rise to tracks in the dipole spectrometer. In 
such cases also the additional tracks are entered 
into the migration matrix thereby taking into ac- 
count the combined effect of the generated parti- 
cle and any secondaries it creates. The absorption 
correction is on average 20%, approximately inde- 
pendent of momentum. Uncertainties in the ab- 
sorption of secondaries in the dipole spectrometer 
material are taken into account by a variation of 
10% of this effect in the simulation. The effect of 
pion decay is treated in the same way as the ab- 
sorption and is 20% at 500 MeV/c and negligible 
at 3 GeV/c. 

The uncertainty in the production of background 
due to tertiary particles is larger. The average cor- 
rection is ps10% and up to 20% at 1 GeV/c. The 
correction includes reinteractions in the detector 
material (mainly carbon) as well as a small com- 
ponent coming from reinteractions in the target. 
The subtraction may be computed by MonteCarlo 
simulations and as most of the encountered mate- 
rial is carbon, the check of the inelastic interac- 
tions of low-energy protons or pions in carbon is 
essential. The validity of the generators used in 
the simulation was checked by the analysis of our 
data with incoming protons and charged pions on 
aluminium and carbon targets at lower momenta 
(3 GeV/c and 5 GeV/c). A 30% uncertainty on 
the secondary production was assumed. 

The average empty-target subtraction amounts to 
«20%. 

Owing to the redundancy of the tracking system 
downstream of the target the detection efficiency 
is very robust under the usual variations of the 
detector performance during the long data tak- 
ing periods. Since the momentum is reconstructed 
without making use of the upstream drift chamber 
module (which is more sensitive in its performance 
to the beam intensity) the reconstruction efficiency 
is uniquely determined by the downstream system. 
No variation of the overall efficiency has been ob- 
served. The performance of the TOFW and CHE 
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system have been monitored to be constant for the 
data taking periods used in this analysis. The cal- 
ibration of the detectors was performed on a day- 
by-day basis. 



B. Error estimation 



The total statistical error takes into account the 
direct error propagation of the statistical errors in 
the raw data and the statistical error which is in- 
curred while obtaining the unfolding matrix from 
the data. The latter component increases the di- 
rect error by a factor two. The procedure is out- 
lined in references [H, [28[. 

Different types of sources induce systematic errors 
for the analysis described here: track yield correc- 
tions (~ 5%), particle identification (~ 0.1%), mo- 
mentum and angular reconstruction (~ 1%) |45j |. 
The strategy to calculate these systematic errors 
and the different methods used for their evaluation 
are described in An additional source of error 
is due to misidentified secondary kaons, which are 
not considered in the particle identification method 
used for this analysis and are subtracted on the 
basis of a Monte Carlo simulation, as in [4[. No 
explicit correction is made for pions coming from 
decays of other particles created in the target, as 
they give a very small contribution according to 
the selection criteria applied in the analysis. As a 
result of these systematic error studies, each error 
source can be represented by a covariance matrix. 
The sum of these matrices describes the total sys- 
tematic error, as explained in Q. 

The experimental uncertainties are shown for a 
typical light target (Be) in Figure for ir + at some 
typical incident beam momenta and in Table HI1 for 
another target (C) at the incident beam momen- 
tum 12 GeV/c. They are very similar for n~ and at 
the other beam energies. Going from lighter (Be, 
C) to heavier targets (Ta, Pb) the corrections for 
7r° (conversion) and absorption/tertiaries increase. 



TABLE II: Summary of the systematic uncertainties 
affecting the computed n + double-differential cross sec- 
tions and the integrated cross-section measurements 
for p-C interactions at 12 GeV/c. The entries of the 
table are weighted bin by bin with the pion production 
yields. 



Error Category 


Sdm (%) 


(%) 


Track yield corrections: 






Reconstruction efficiency 


1.1 


0.5 


Pion, proton absorption 


3.7 


3.2 


Tertiary subtraction 


8.6 


3.7 


Empty target subtraction 


1.2 


1.2 


Sub-total 


9.5 


5.1 


Particle Identification: 






Electron veto 


<0.1 


<0.1 


Pion, proton ID correction 


0.1 


0.1 


Kaon subtraction 


<0.1 


<0.1 


Sub-total 


0.1 


0.1 


Momentum reconstruction: 






Momentum scale 


2.8 


0.3 


Momentum resolution 


0.8 


0.3 


Sub-total 


2.9 


0.4 


Angle reconstruction: 






Angular scale 


1.3 


0.5 


Total syst. 


10.0 


5.1 


Overall normalization: 


2.0 


2.0 



same order. 



The overall normalization has an uncertainty of 
~ 2% and is mainly due to the uncertainty in 
the efficiency that beam protons counted in the 
normalization actually hit the target, with smaller 
components from the target density and the beam 
particle counting procedure. 

On average the total integrated systematic error 
is around 5 — 6%, with a differential bin to bin 
systematic error of the order of 10 — 11%, to be 
compared with a statistical integrated (bin-to-bin 
differential) error of - 2 - 3% (~ 10 - 13%). Sys- 
tematic and statistical errors are roughly of the 
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tion+tertiares interactions, black dotted line for track efficiency and target pointing efficiency (mostly < 1%), 
black solid line for momentum scale+resolution and angle scale for a typical target (Be) at some typical incident 
momenta. 
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IV. EXPERIMENTAL RESULTS 



The measured double-differential cross-sections for 
the production of tt + and 7r~ in the laboratory sys- 
tem as a function of the momentum and the po- 
lar angle for each incident beam momentum are 
shown in Figures [3] to [9] for solid targets from Be 
to Pb. The error bars shown are the square-roots 
of the diagonal elements in the covariance matrix, 
where statistical and systematic uncertainties are 
combined in quadrature. The correlation of the 
statistical errors (introduced by the unfolding pro- 
cedure) are typically smaller than 20% for adja- 
cent momentum bins and even smaller for adjacent 
angular bins. The correlations of the systematic 
errors are larger, typically 80% for adjacent bins. 
The overall normalization error (< 2%) is not in- 
cluded in the error bars j4^|. The results of this 
analysis are also fully tabulated and reported in 
appendix [SJ 

An overall fit with a Sanford-Wang parametriza- 
tion [29| has been done by using all solid targets 
and available beam momenta data (p > 5 GeV/c), 
see Section llVAl for details, and is shown as a solid 
line on all figures. 

The dependence of the averaged pion yields on 
the incident beam momentum is shown in Fig. 1101 
The 7r~ yields, averaged over two angular re- 
gions (0.05 rad < 9 < 0.15 rad and 0.15 rad < 
9 < 0.25 rad) and four momentum regions 
(0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 
2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 
3.5 GeV/c < p < 4.5 GeV/c), are shown in the 
left panel and the tt + data averaged over the same 
regions in the right panel, for four different beam 
momenta. Whereas the beam energy dependence 
of the yields in the p-Be, p-C data differs clearly 
from the dependence in the p-Ta, p-Pb data one 
can observe that the p-Al, p-Cu and p-Sn data 
display a smooth transition between them. The 
dependence in the p-Be, p-C data is much more 
flat with a saturation of the yield between 8 GeV/c 
and 12 GeV/c with the p-Al, p-Cu and p-Sn show- 
ing an intermediate behaviour. 

The dependence of the averaged pion yields on 
the atomic number A is instead shown in Fig. 111! 
The 7T~ yields, averaged over two angular re- 
gions (0.05 rad < 9 < 0.15 rad and 0.15 rad < 
9 < 0.25 rad) and four momentum regions 
(0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 
2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 
3.5 GeV/c < p < 4.5 GeV/c), are shown in the 
left panel and the ir + data averaged over the same 
regions in the right panel, for four different beam 
momenta. One observes a smooth behaviour of 
the averaged yields. The A-dependence is slightly 



different for 7r~ and ir + production, the latter satu- 
rating earlier towards higher A, especially at lower 
beam momenta. 
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FIG. 3: p-Be differential cross sections: left panel n~ production, right panel n + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 4: p-C differential cross sections: left panel ir~ production, right panel n + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 5: p-Al differential cross sections: left panel 7r~ production, right panel 7r + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 6: p-Cu differential cross sections: left panel ir~ production, right panel 7r + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 7: p-Sn differential cross sections: left panel 7r~ production, right panel 7r + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
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FIG. 8: p-Ta differential cross sections: left panel 7r~ production, right panel ir + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 9: p-Pb differential cross sections: left panel 7r~ production, right panel n + production. The curves represent 
the global parametrization as described in the text. In the top right corner of each plot the covered angular range 
is shown in mrad. 
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FIG. 10: The dependence on the beam momentum of the 7r (top panel) and tt + (bottom panel) production 
yields in p-Be, p-C, p-Al, p-Cu, p-Sn, p-Ta, p-Pb interactions averaged over two forward angular regions 
(0.05 rad < 9 < 0.15 rad and 0.15 rad < 9 < 0.25 rad) and four momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 
1.5 GeV/c < p < 2.5 GeV/c, 2.5 GeV/c < p < 3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for the four 
different incoming beam energies. Data points for different target nuclei and equal momenta are slightly shifted 
horizontally with respect to each other to increase the visibility. 
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FIG. 11: The dependence on the atomic number A of the pion production yields in p-Be, p-Al, p-C, p- 
Cu, p-Sn, p-Ta, p-Pb interactions averaged over two forward angular regions (0.05 rad < 6 < 0.15 rad and 
0.15 rad < 6 < 0.25 rad) and four momentum regions (0.5 GeV/c < p < 1.5 GeV/c, 1.5 GeV/c < p < 2.5 GeV/c, 
2.5 GeV/c < p < 3.5 GeV/c and 3.5 GeV/c < p < 4.5 GeV/c), for the four different incoming beam energies. 
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A. Pion production data parametrization 



At low energies, it is common to use the empir- 
ical data parametrization for pion production in 
proton-nucleus interactions originally developed by 
Sanford and Wang 29]. This parametrization has 
the functional form: 



d 2 a(pA 



r±X) 



dpdfl 



ciexp[B]p C2 (l 



Pbe 



(2) 



where: 



B = 



~ C 3 — 
Pbe: 



c 6 0(p - c 7 p bcam cos C8 6) (3) 



X denotes any system of other particles in the 
final state, Pbeam is the proton beam momen- 
tum in GeV/c, p and 8 are the momentum 
and angle in units of GeV/c and radians, re- 
spectively, d 2 a/(dpdil) is expressed in units of 
mb/(GeV/c)/sr and the parameters ci, . . . , c% are 
obtained from fits to the meson production data. 

The parameter c\ is an overall normalization fac- 
tor, the four parameters 02,03,04,05 describe the 
momentum distribution of the secondary pions in 
the forward direction, and the three parameters 
06,07,08 describe the corrections to the pion mo- 
mentum distribution for pion production angles 
that are different from zero. The ir^ production 
data on solid targets reported here have been fit- 
ted simultaneously to the empirical Sanford- Wang 
formula. In the % 2 minimization, the full error ma- 
trix was used. 

To go from a baseline nuclear targets (typically Be) 
to another nuclear target A, a correction factor 



corr = (A/A Be ) a 



(4) 



was introduced, with a — + a\ x xf + 012 x x F 
where xp is the Feynman x, see references [loj and 
[30j for details. A 9-parameter (11-parameter) fit 
was done over 24 (25) tt~(it + ) datasets [13], corre- 
sponding to 1440 (1472) experimental points. The 
goodness-of-fit of the Sanford- Wang parametriza- 
tion hypothesis for the HARP results can be as- 
sessed by considering the best-fit \ 2 value of 
x 2 m . m = 13030 (8061) for 1431 (1461) degrees of 
freedom for the 7r _ (tt + ) production, indicating 
a very poor fit quality. In particular, inspection 
of the HARP inclusive pion production double- 
differential cross-section, and resulting Sanford- 
Wang parametrization, points to a description of 
the ratio g{0) of the pion momentum distribution 
at 6 ^ with respect to the 9 = pion momen- 
tum distribution that is more complicated than 



TABLE III: Sanford- Wang parameters and errors ob- 
tained by fitting the n + (tv~) datasets. The errors 
refer to the 68.27% confidence level for eleven pa- 
rameters (Ax 2 = 12.6) for 7r + and nine parameters 
(Ax 2 = 10.4) I for n~. 



Parameter 



Cl 


(381.3 ±40.5) 


(307.6 ± 19.4) 


C2 


(0.88 ± 0.07) 


(0.57 ± 0.06) 


C3 


(9.16 ±0.95) 


(27.36 ± 1.06) 


C4 


(1.38 ±0.09) 


(1.86 ±0.04) 


C5 


(1.66 ±0.12) 


(2.23 ± 0.04) 


C6 


(3.62 ±0.14) 


(3.04 ± 0.08) 


C7 


(0.05 ± 0.04) 




C8 


(128.6 ±61.8) 




Q() 


(0.69 ± 0.04) 


(0.72 ± 0.04) 


ai 


(-0.91 ±0.21) 


(-1.36 ±0.20) 


a 2 


(0.34 ±0.21) 


(2.18 ±0.21) 



what can be accommodated within the Sanford- 
Wang formula, where this ratio is given by g(0) = 
exp[-c 6 6»(]j - p c )], with p c = c 7 p boam cos C8 6. 

The overall fit may be used as a fast approxima- 
tion of HARP data valid within a factor 2-3 of the 
quoted experimental errors. The best-fit values of 
the parameters are reported in Table ITTT1 together 
with their errors. The fit parameter errors are esti- 



mated by requiring A\ 



X 



-X 2 ■ 



12.6 (10.4), 



corresponding to the 68.27% confidence level re- 
gion for eleven (nine) variable parameters. Signif- 
icant correlations among fit parameters are found, 
as shown by the correlation matrix given in Ta- 
bles ED and El 

To show the trend of the Sanford- Wang global fit 
of all HARP datasets, figure [T^] reports the com- 
parison, at 8 GeV/c and 12 GeV/c, between pion 
production data and the above parametrization. 

For the 8.9 GeV/c MiniBooNE/SciBooNE beam- 
line and the 12.9 GeV/c K2K beamline two ad- 
hoc Sanford- Wang parametrizations, using only 
the relevant HARP datasets, have been published 
in references [l[ and • Given the poor description 
of HARP pion production data in terms of the orig- 
inal Sanford- Wang parametrization, one extra pa- 
rameter to better describe the angular dependence 
was introduced in the fit reported in reference @ ■ 
In the global fit presented here, this extra param- 
eter is not used for simplicity, as it was found that 
it did not improve the fit quality in a significant 
way. 

As a final remark, we stress again that due to the 
poor fit quality our global fit may be just consid- 
ered as a simple way to summarize an extended set 



TABLE IV: Correlation coefficients among the Sanford-Wang fit parameters, obtained by fitting the data for 
production. 



Parameter ci C2 C3 C4 C5 ce c-j c% ao ai oti 

ci 1.000 

c 2 0.388 1.000 

c 3 0.089 -0.349 1.000 

c 4 -0.725 -0.485 0.103 1.000 

c 5 -0.580 -0.651 0.585 0.840 1.000 

c 6 0.309 0.358 0.069 -0.262 -0.106 1.000 

c 7 -0.612 -0.453 -0.237 0.299 0.137 -0.532 1.000 

cs -0.017 0.187 0.027 -0.080 -0.017 0.513 -0.319 1.000 

a -0.528 0.131 -0.030 0.364 0.289 -0.008 0.013 -0.025 1.000 

ai 0.451 -0.095 -0.061 -0.529 -0.429 0.065 -0.017 0.035 -0.827 1.000 

q 2 -0.377 0.059 0.228 0.557 0.498 -0.069 0.012 -0.027 0.672 -0.945 1.000 



TABLE V: Correlation coefficients among the Sanford-Wang fit parameters, obtained by fitting the data for 
production. 



Parameter ci C2 C3 C4 C5 cq ao ai 012 

ci 1.000 

c 2 -0.151 1.000 

c 3 0.040 -0.517 1.000 

c 4 -0.544 -0.263 0.047 1.000 

c 5 -0.441 -0.567 0.625 0.794 1.000 

c 6 0.134 0.181 0.365 -0.269 0.027 1.000 

a -0.728 0.307 -0.049 0.366 0.266 0.010 1.000 

qi 0.492 -0.325 0.037 -0.493 -0.318 0.069 -0.809 1.000 

q 2 -0.342 0.285 0.040 0.492 0.330 -0.073 0.636 -0.950 1.000 



of data (~ 1000 experimental data points) using a 
formula with about ten parameters. This may be 
useful in the intial phase of an experiment design 
or a Monte Carlo validation. 
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FIG. 12: Comparison of HARP production data at 8 GeV/c and 12 GeV/c with the Sanford-Wang global fit. 
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B. Comparison with Monte Carlo generators. 

In the following we will show only some compar- 
isons with two widely available Monte Carlo gen- 
erators: GEANT4 version 7.1 [27[ and MARS ver- 
sion 15.07 [3l[, using different models. The com- 
parison will be shown for a limited set of plots and 
only for the Be and Ta targets, as examples of a 
light and a heavy target. In both generators, no 
single model is applicable to all energies and a tran- 
sition between low energy models and high energy 
models, at about 5 GeV/c-10 GeV/c, has to be 
done. 

At intermediate energies (up to 5 GeV/c- 
10 GeV/c), GEANT4 uses two types of intra- 
nuclear cascade models: the Bertini model [32, 33] 
(valid up to ~10 GeV/c) and the Binary model [34j 
(valid up to ^3 GeV/c). Both models treat the tar- 
get nucleus in detail, taking into account density 
variations and tracking in the nuclear field. The 
Binary model is based on hadron collisions with 
nucleons, giving resonances that decay according 
to their quantum numbers. The Bertini model is 
based on the cascade code reported in [35] and 
hadron collisions are assumed to proceed accord- 
ing to free-space partial cross sections corrected 
for nuclear field effects and final state distributions 
measured for the incident particle types. 

At higher energies, instead, two parton string mod- 
els, the quark-gluon string (QGS) model [H, |3(| 
and the Fritiof (FTP) model (36| are used, in addi- 
tion to a High Energy Parametrized model (HEP) 
derived from the high energy part of the GHEISHA 
code used inside GEANT3 (37| . The parametrized 
models of GEANT4 (HEP and LEP) are intended 
to be fast, but conserve energy and momentum on 
average and not event by event. 

A realistic GEANT4 simulation is built by com- 
bining models and physics processes into what is 
called a "physics list" . In high energy calorime- 
try the two most commonly used are the QGSP 
physics list, based on the QGS model, the pre- 
compound nucleus model and some of the Low 
Energy Parametrized (LEP) model ,48] and the 
LHEP physics list [H| based on the parametrized 
LEP and HEP models. 

The MARS code system [3l[ uses instead as 
basic model an inclusive approach multiparticle 
production originated by R. Feynmann. Above 
3 GeV phenomenological particle production mod- 
els are used, while below 5 GeV a cascade-exciton 
model (39l ] combined with the Fermi break-up 
model, the coalescence model, an evaporation 
model and a multifragmentation extension are 
used instead. 



In figures [T31 to [20] data have been compared to 
Monte Carlo predictions, using beryllium and tan- 
talum as examples of a light and a heavy target. 

Computing the \ 2 between models and data them- 
selves, where a systematic error going from 0% to 
50% have been added to simulation results we ob- 
tain the results shown in Table IVIl Normaliza- 
tion factors, between data and Monte Carlo simu- 
lations, are shown instead in Table rVlTl 

Over the full energy range covered by the HARP 
experiment, the best comparison is obtained with 
the MARS Monte Carlo. This may be owing to the 
fact that MARS is using different models in differ- 
ent energy regions, equivalent to use a collection 
of models as implemented in the "physics lists" of 
GEANT4. 

The full set of HARP data, taken with targets 
spanning the full periodic table of elements, with 
small total errors and full coverage of the solid an- 
gle in a single detector may help the validation of 
models used in hadronic simulations in the difficult 
energy range between 3 and 15 GeV/c of incident 
momentum, as done in reference [40| . 
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FIG. 13: Comparison of HARP double-differential ir^ cross sections for p-Be at 3 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): binary model grey line, Bertini 
model black solid line, LHEP model dotted line, MARS model grey solid line. 
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FIG. 14: Comparison of HARP double-differential ^ cross sections for p-Be at 5 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): Binary model grey line, Bertini 
model black solid line, LHEP model dotted line, MARS model grey solid line. 
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TABLE VI: Computed \ 2 between data and MonteCarlo simulations, assuming a 0%(50%) systematica on simu- 
lation 





3 GeV 


5 GeV 


8 GeV 


12 GeV 






TV 










7V+ 


TT 


ndof 


16 


16 


32 


32 


36 


36 


40 


40 


model 


Beryllium 


Bcrtini 


725.0 (39.6) 


602.7 (35.6) 


2334.5 (74.0) 


1722.0 (79.9) 


3372.7 (171.0) 


2461.4 (68.6) 


3501.3 (331.9) 


2390.8 (90.3) 


Binary 


690.1 (39.7) 


1011.1 (41.6) 


1020.8 (65.6) 


3708.8 (91.8) 


1411.6 (74.0) 


3481.4 (77.3) 






LHEP 


428.4 (36.8) 


103.7 (14.9) 


265.4(34.8) 


1109.7 (71.8) 


539.5 (29.1) 


558.2 (31.4) 


893.2 (34.6) 


1084.4 (33.4) 


QGSP 










226.6 (35.0) 


740.3 (49.4) 


809.7 (31.3) 


972.8 (33.3) 


FTFP 












1805.7 (65.7) 


196.2 (16.6) 


1113.6 (36.3) 


MARS 


12.3 (5.1) 


91.7 (17.2) 


33.9 (10.1) 


17.0 (9.1) 


65.4 (9.0) 


122.3 (16.0) 


62.6 (7.5) 


70.0 (7.4) 


model 


Tantalum 


Bcrtini 


493.1 (30.6) 


380.6 (25.9) 


517.5 (41.3) 


2144.2 (80.2) 


1186.1 (84.5) 


735.6 (45.1) 


1455.8 (162.4) 


709.7 (68.4) 


Binary 


901.3 (35.2) 


784.9 (30.3) 


743.4 (46.3) 


4075.1 (82.4) 


1093.1 (53.0) 


2535.0 (53.8) 






LHEP 


600.7 (28.7) 


57.5 (7.7) 


201.0 (28.8) 


1022.6 (61.7) 


371.4 (22.3) 


516.3 (21.7) 


434.7 (14.8) 


884.9 (25.4) 


QGSP 










439.9 (43.2) 


678.1 (36.9) 


213.3 (12.5) 


636.0 (18.5) 


FTFP 












365.7 (23.6) 


135.5 (31.2) 


259.9 (15.7) 


MARS 


145.8 (14.3) 


70.5 (9.1) 


143.2 (32.1) 


94.4 (29.5) 


121.1 (18.1) 


53.3 (6.7) 


101.2 (7.6) 


54.0 (5.8) 



TABLE VII: Normalization factors data-simulation. 



model 


Be 3 GeV 


Ta 3 GeV 


Be 5 GeV 


Ta 5 GeV 


Be 8 


GeV 


Ta 8 GeV 


Be 12 GeV 


Ta 12 GeV| 




TT + 


TT 


TT+ 


TT 


TT + 


TT 


TT + 


TT 


tt+ 




TT + 


TT 


TT + 


TT 


TT + 


TT 


Bertini 


0.35 


1.02 


0.45 


0.53 


0.70 


1.12 


0.29 


0.35 


1.22 


1.54 


0.84 


1.08 


1.75 


1.81 


1.27 


1.50 


Binary 


0.36 


0.75 


0.28 


0.34 


0.73 


0.88 


0.16 


0.23 


0.99 


1.05 


0.50 


0.56 










LHEP 


0.40 


0.86 


0.81 


0.91 


0.76 


0.98 


0.36 


0.45 


0.78 


0.91 


0.58 


0.66 


0.75 


0.82 


0.54 


0.59 


QGSP 


















1.40 


1.43 


0.80 


0.75 


0.80 


0.88 


0.64 


0.67 


FTFP 






















0.46 


0.65 


1.00 


1.10 


0.63 


0.77 


MARS 


0.83 


1.29 


1.10 


1.16 


1.21 


1.38 


1.17 


1.35 


1.10 


1.21 


0.90 


0.85 


1.02 


1.02 


0.92 


0.82 
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FIG. 15: Comparison of HARP double-differential cross sections for p-Be at 8 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): QGSC model dotted line, QGSP 
model black solid line, MARS model dashed line. 
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FIG. 16: Comparison of HARP double-differential % ± cross sections for p-Be at 12 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details) :FTFP model dotted line, QGSP 
model black solid line, MARS dashed line 
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FIG. 17: Comparison of HARP double-differential ^ cross sections for p-Ta at 3 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): Binary model grey line, Bertini 
model black solid line, LHEP model dotted line, MARS model grey solid line. 
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FIG. 18: Comparison of HARP double-differential ty cross sections for p-Ta at 5 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): Binary model grey line, Bertini 
model black solid line, LHEP model dotted line, MARS model grey solid line. 
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FIG. 19: Comparison of HARP double-differential ^ cross sections for p-Ta at 8 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details): QGSC model dotted line, QGSP 
model black solid line, MARS dashed line. 
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FIG. 20: Comparison of HARP double-differential cross sections for p-Ta at 12 GeV/c with GEANT4 and 
MARS MC predictions, using several generator models (see text for details) : FTFP model dotted grey line, QGSP 
model black solid line, MARS dashed line. 
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V. SUMMARY AND CONCLUSIONS 



In this paper we report our final results on double- 
differential cross sections for the production of pos- 
itive and negative pions in the kinematic range 0.5 
GeV/c< p„ < 8 GeV/c and 0.025 rad < 6 n < 0.25 
rad from the collisions of protons of 3, 5, 8 and 12 
GeV/c on beryllium, aluminium, carbon, copper, 
tin, tantalum targets of 5% A/. In addition results 
at 8.9 GeV/c (Be only), 12.9 GeV/c (Al only) and 
12 GeV/c (N2 and O2) are reported. 

A parametrization, inspired by the Sanford-Wang 
formula, of all our datasets at all energies is also 
presented. This may be used as a fast approxima- 
tion to HARP data, valid within a factor 2-3 of the 
quoted errors. 

The pion yield averaged over different momentum 
and angular ranges increase smoothly with the 
atomic number A of the target and with the energy 
of the incoming proton beam. The A-dependence 
is slightly different for it~ and ir + production, the 
latter saturating earlier towards higher A, espe- 
cially at lower beam momenta. 

Comparisons with GEANT4 and MARS genera- 
tors are presented. 
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APPENDIX A: CROSS-SECTION DATA 



The following tables report the measured differ- 
ential cross-section for positive and negative pion 
production in interactions of 3, 5, 8 and 12 GeV/c 
momentum protons on different types of nuclear 
targets. The data are presented in the kinematic 
range of 0.5 GeV/c < p n < 8 GeV/c and 0.05 rad 
< n < 0.25 rad. The overall normalization un- 
certainty (< 2%) is not included in the reported 
errors. 

Results at higher incoming beam momenta, from 
8 GeV/c to 12 GeV/c, are also presented extending 
to a lower value of the polar angle 9 (0.025 rad), 
using a finer binning. 

The cross-section values in some bins of the follow- 
ing tables are omitted and replaced by the symbol 
* ± *, due to instabilities in the unfolding proce- 
dure coming from the low statistics available. This 
appears especially evident for the 3 GeV/c data 
where the total negative final state pion statistics 
is often well below one hundred events. 
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TABLE VIII: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 er 7r / (dpdQ) , for p-Be interactions at 3, 5, 8, 8.9, 12 GeV/c. Each row refers to a different (p mm < p < 
Pmax,# mm < < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central 
value as well as the square-root of the diagonal elements of the covariance matrix are given. 
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0.158±0 


011 






3.50 


4.00 


± 


0.003±0 


002 


0.028±0.005 


0.065±0.006 


0.122±0 


010 






4.00 


5.00 


± 


0.001±0 


001 


0.021±0.004 


0.036±0.004 


0.073±0 


006 






5.00 


6.50 


± 


± 




0.003±0.001 


0.009±0.002 


0.033±0 


003 






6.50 


8.00 


± 


± 




± 


± 


0.011±0 


002 


0.100 


0.150 


0.50 


1.00 


0.10±0.03 


0.13±0 


02 


0.18±0.02 


0.17±0.02 


0.20±0 


02 






1.00 


1.50 


0.028±0.013 


0.13±0 


02 


0.17±0.02 


0.19±0.02 


0.22±0 


02 






1.50 


2.00 


O.OOliO.001 


0.066±0 


011 


0.181±0.015 


0.168±0.013 


0.22±0 


02 






2.00 


2.50 




0.035±0 


007 


0.092±0.010 


0.130±0.010 


0.16±0 


02 






2.50 


3.00 


± 


0.011±0 


005 


0.049±0.006 


0.075±0.007 


0.107±0 


010 






3.00 


3.50 


± 


0.001±0 


001 


0.023±0.004 


0.040±0.005 


0.078±0 


008 






3.50 


4.00 


± 






0.012±0.002 


0.024±0.003 


0.055±0 


006 






4.00 


5.00 


± 






0.004±0.001 


0.010±0.002 


0.033±0 


004 






5.00 


6.50 


± 


± 




*_j_=f= 




0.006±0 


001 






6.50 


8.00 


± 


± 




± 


± 






0.150 


0.200 


0.50 


1.00 


0.07±0.02 


0.10±0 


02 


0.14±0.02 


0.17±0.02 


0.16±0 


02 






1.00 


1.50 


0.010±0.007 


0.11±0 


02 


0.143±0.013 


0.155±0.012 


0.19±0 


02 






1.50 


2.00 


*_|_* 


0.029±0 


007 


0.080±0.008 


0.102±0.008 


0.129±0 


013 






2.00 


2.50 


*_|_* 


0.013±0 


005 


0.055±0.007 


0.065±0.007 


0.091±0 


011 






2.50 


3.00 


± 


0.002±0 


002 


0.020±0.005 


0.034±0.004 


0.047±0 


007 






3.00 


3.50 


± 


0.001±0 


001 


0.008±0.002 


0.020±0.002 


0.024±0 


004 






3.50 


4.00 


± 


*_j_* 




0.003±0.001 


0.011±0.002 


0.014±0 


003 






4.00 


5.00 


± 






O.OOliO.001 


0.005±0.001 


0.008±0 


002 






5.00 


6.50 


± 


± 






*_|_=f= 


0.002±0 


001 






6.50 


8.00 


± 


± 




± 


± 






0.200 


0.250 


0.50 


1.00 


0.05±0.02 


0.09±0 


02 


0.14±0.02 


0.120±0.015 


0.12±0 


02 






1.00 


1.50 


0.017±0.011 


0.032±0 


009 


0.082±0.012 


0.071±0.010 


0.080±0 


013 






1.50 


2.00 


=f=_|_=f= 


0.027±0 


009 


0.057±0.010 


0.049±0.007 


0.064±0 


011 






2.00 


2.50 


*_|_* 


0.016±0 


006 


0.038±0.007 


0.033±0.006 


0.053±0 


010 






2.50 


3.00 


± 


0.003±0 


002 


0.025±0.005 


0.025±0.003 


0.031±0 


007 






3.00 


3.50 


± 


0.002±0 


002 


0.011±0.003 


0.014±0.002 


0.015±0 


004 






3.50 


4.00 


± 


0.001±0 


002 


0.005±0.002 


0.006±0.001 


0.007±0 


002 






4.00 


5.00 


± 


0.001±0 


001 


0.003±0.001 


0.004±0.001 


0.005±0 


002 






5.00 


6.50 


± 


± 




O.OOliO.001 




0.001±0 


001 






6.50 


8.00 


± 


± 




± 


± 


0.000±0 


002 
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TABLE IX: HARP results for the double-differential 7r~ production cross-section in the laboratory system, 
d 2 er 7r / (dpdQ,) , for p-Be interactions at 3, 5, 8, 8.9, 12 GeV/c. Each row refers to a different (p mm < p < 
Pma.x,8 m [ n < < #max) bin, where p and 8 are the pion momentum and polar angle, respectively. The central 
value as well as the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


JJmax 




d 2 a n /(dpdQ,) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


8.9 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.017±0.011 


0.07±0.02 


0.12±0.02 


0.118±0.014 


0.14±0.02 






1.00 


1.50 


0.002±0.003 


0.052±0.010 


0.116±0.011 


0.141±0.010 


0.158±0.014 






1.50 


2.00 


0.001±0.002 


0.030±0.007 


0.101±0.009 


0.118±0.008 


0.159±0.013 






2.00 


2.50 




0.023±0.006 


0.060±0.006 


0.082±0.006 


0.124±0.011 






2.50 


3.00 


± 


0.011±0.004 


0.033±0.006 


0.049±0.005 


0.097±0.009 






3.00 


3.50 


± 


0.001±0.001 


0.019±0.003 


0.037±0.003 


0.076±0.007 






3.50 


4.00 


± 




0.012±0.002 


0.024±0.004 


0.045±0.005 






4.00 


5.00 


± 




0.007±0.002 


0.013±0.002 


0.034±0.004 






5.00 


6.50 


± 


± 




0.002±0.001 


0.010±0.002 






6.50 


8.00 


± 


± 


± 


± 


0.003±0.001 


0.100 


0.150 


0.50 


1.00 


0.04±0.02 


0.09±0.02 


0.15±0.02 


0.15±0.02 


0.19±0.02 






1.00 


1.50 


0.007±0.006 


0.054±0.011 


0.093±0.010 


0.126±0.011 


0.152±0.015 






1.50 


2.00 


0.001±0.002 


0.022±0.006 


0.085±0.008 


0.085±0.008 


0.115±0.012 






2.00 


2.50 


*_|_* 


0.008±0.003 


0.050±0.006 


0.056±0.006 


0.085±0.009 






2.50 


3.00 


± 


0.006±0.003 


0.024±0.004 


0.036±0.004 


0.062±0.008 






3.00 


3.50 


± 


0.001±0.001 


0.014±0.002 


0.020±0.002 


0.048±0.006 






3.50 


4.00 


± 




0.008±0.002 


0.015±0.002 


0.027±0.004 






4.00 


5.00 


± 




0.003±0.001 


0.005±0.001 


0.014±0.002 






5.00 


6.50 


± 


± 




0.001±0.000 


0.004±0.001 






6.50 


8.00 


± 


± 


± 


± 




0.150 


0.200 


0.50 


1.00 


0.012±0.008 


0.09±0.02 


0.13±0.02 


0.16±0.02 


0.18±0.02 






1.00 


1.50 


0.006±0.005 


0.044±0.010 


0.086±0.010 


0.093±0.009 


0.110±0.013 






1.50 


2.00 


0.007±0.007 


0.017±0.005 


0.048±0.007 


0.060±0.007 


0.083±0.010 






2.00 


2.50 




0.008±0.003 


0.037±0.005 


0.038±0.005 


0.064±0.008 






2.50 


3.00 


± 




0.018±0.003 


0.018±0.003 


0.033±0.005 






3.00 


3.50 


± 




0.009±0.002 


0.010±0.002 


0.025±0.004 






3.50 


4.00 


± 




0.002±0.001 


0.005±0.001 


0.016±0.003 






4.00 


5.00 


± 






0.002±0.001 


0.007±0.002 






5.00 


6.50 


± 


± 




*_|_=f= 








6.50 


8.00 


± 


± 


± 


± 




0.200 


0.250 


0.50 


1.00 


0.03±0.02 


0.07±0.02 


0.096±0.014 


0.102±0.014 


0.13±0.02 






1.00 


1.50 


0.001±0.003 


0.030±0.009 


0.081±0.011 


0.090±0.011 


0.106±0.015 






1.50 


2.00 


*_|_=f= 


0.017±0.008 


0.064±0.012 


0.077±0.012 


0.11±0.02 






2.00 


2.50 


*_|_* 


0.004±0.003 


0.029±0.006 


0.040±0.007 


0.061±0.011 






2.50 


3.00 


± 


0.000±0.001 


0.008±0.003 


0.019±0.004 


0.028±0.007 






3.00 


3.50 


± 




0.003±0.002 


0.008±0.002 


0.013±0.004 






3.50 


4.00 


± 




0.002±0.001 


0.003±0.001 


0.006±0.002 






4.00 


5.00 


± 




O.OOliO.001 


0.001±0.000 


O.OOliO.001 






5.00 


6.50 


± 


± 


O.OOOiO.001 


*_|_* 








6.50 


8.00 


± 


± 


± 


± 


O.OOOiO.001 
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TABLE X: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 cr 7r / (dpdfl) , for p-C interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < V < Pmax,# mm < 
9 < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


^min 


Pmax 




d 2 a' K+ /(dpdQ) 




(rad) 


(rad) (GeV/c) 


(GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.11±0.03 


0.12±0.02 


0.18±0.02 


0.19±0.03 






1.00 


1.50 


0.015±0.012 


0.18±0.02 


0.20±0.02 


0.24±0.02 






1.50 


2.00 


0.013±0.007 


0.113±0.014 


0.25±0.02 


0.28±0.02 






2.00 


2.50 


O.OOliO.001 


0.044±0.008 


0.185±0.014 


0.27±0.02 






2.50 


3.00 


± 


0.023±0.005 


0.135±0.013 


0.23±0.02 






3.00 


3.50 


± 


0.013±0.005 


0.060±0.010 


0.167±0.014 






3.50 


4.00 


± 


0.003±0.002 


0.030±0.007 


0.145±0.014 






4.00 


5.00 


± 




0.021±0.004 


0.084±0.007 






5.00 


6.50 


± 


± 


0.004±0.001 


0.037±0.004 






6.50 


8.00 


± 


± 


± 


0.016±0.003 


0.100 


0.150 


0.50 


1.00 


0.12±0.04 


0.16±0.02 


0.21±0.02 


0.21±0.03 






1.00 


1.50 


0.024±0.015 


0.15±0.02 


0.22±0.02 


0.23±0.02 






1.50 


2.00 




0.094±0.013 


0.20±0.02 


0.26±0.03 






2.00 


2.50 




0.022±0.005 


0.123±0.014 


0.21±0.02 






2.50 


3.00 


± 


0.012±0.004 


0.049±0.007 


0.135±0.014 






3.00 


3.50 


± 


0.002±0.001 


0.027±0.006 


0.093±0.011 






3.50 


4.00 


± 


*_|_* 


0.012±0.003 


0.061±0.008 






4.00 


5.00 


± 


*_|_* 


0.004±0.001 


0.037±0.005 






5.00 


6.50 


± 


± 




0.010±0.002 






6.50 


8.00 


± 


± 


± 


0.003±0.001 


0.150 


0.200 


0.50 


1.00 


0.09±0.03 


0.16±0.02 


0.22±0.03 


0.24±0.03 






1.00 


1.50 


0.019±0.012 


0.094±0.014 


0.19±0.02 


0.21±0.02 






1.50 


2.00 




0.042±0.008 


0.108±0.011 


0.15±0.02 






2.00 


2.50 




0.019±0.006 


0.060±0.009 


0.091±0.012 






2.50 


3.00 


± 


0.005±0.003 


0.032±0.005 


0.050±0.008 






3.00 


3.50 


± 


=f=_|_* 


0.019±0.004 


0.031±0.006 






3.50 


4.00 


± 


*_|_* 


0.007±0.002 


0.024±0.005 






4.00 


5.00 


± 


*_|_* 


0.002±0.001 


0.011±0.003 






5.00 


6.50 


± 


± 




0.004±0.001 






6.50 


8.00 


± 


± 


± 


*_|_=f= 


0.200 


0.250 


0.50 


1.00 


0.07±0.03 


0.12±0.02 


0.17±0.02 


0.17±0.03 






1.00 


1.50 


0.012±0.021 


0.053±0.012 


0.12±0.02 


0.08±0.02 






1.50 


2.00 


O.OOliO.010 


0.024±0.007 


0.076±0.013 


0.09±0.02 






2.00 


2.50 




0.010±0.004 


0.055±0.010 


0.044±0.010 






2.50 


3.00 


± 


0.007±0.003 


0.032±0.007 


0.028±0.008 






3.00 


3.50 


± 


0.002±0.001 


0.017±0.005 


0.018±0.005 






3.50 


4.00 


± 


0.000±0.001 


0.007±0.003 


0.008±0.004 






4.00 


5.00 


± 


*_|_=f= 


0.003±0.002 


0.005±0.003 






5.00 


6.50 


± 


± 


0.001±0.001 


0.002±0.001 






6.50 


8.00 


± 


± 


± 


0.001±0.003 
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TABLE XI: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 er 7r / (dpdfl) , for p-C interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < P < Pmax,# mm < 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


^min 


Pmax 




d 2 a n /{dpdQ) 




(rad) 


(rad) (GeV/c) 


(GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.02±0.02 


0.07±0.02 


0.13±0.02 


0.13±0.02 






1.00 


1.50 


0.006±0.008 


0.062±0.010 


0.16±0.02 


0.20±0.02 






1.50 


2.00 


0.002±0.003 


0.038±0.008 


0.108±0.011 


0.19±0.02 






2.00 


2.50 




0.015±0.004 


0.071±0.008 


0.148±0.015 






2.50 


3.00 


± 


0.005±0.002 


0.032±0.007 


0.102±0.013 






3.00 


3.50 


± 


*_|_* 


0.024±0.005 


0.077±0.009 






3.50 


4.00 


± 


*_|_* 


0.025±0.005 


0.056±0.008 






4.00 


5.00 


± 


*_|_* 


0.008±0.002 


0.041±0.005 






5.00 


6.50 


± 


± 




0.009±0.002 






6.50 


8.00 


± 


± 


± 


0.002±0.001 


0.100 


0.150 


0.50 


1.00 


0.01±0.01 


0.10±0.02 


0.20±0.03 


0.24±0.03 






1.00 


1.50 


0.002±0.004 


0.047±0.010 


0.144±0.015 


0.21±0.02 






1.50 


2.00 




0.019±0.005 


0.080±0.009 


0.15±0.02 






2.00 


2.50 




0.013±0.004 


0.050±0.007 


0.095±0.012 






2.50 


3.00 


± 


0.006±0.002 


0.028±0.005 


0.079±0.010 






3.00 


3.50 


± 


0.001±0.001 


0.015±0.003 


0.050±0.008 






3.50 


4.00 


± 


*_|_* 


0.009±0.002 


0.028±0.004 






4.00 


5.00 


± 


*_|_* 


0.003±0.001 


0.017±0.004 






5.00 


6.50 


± 


± 




0.003±0.001 






6.50 


8.00 


± 


± 


± 




0.150 


0.200 


0.50 


1.00 


0.03±0.02 


0.079±0.015 


0.16±0.02 


0.22±0.03 






1.00 


1.50 


0.006±0.007 


0.041±0.009 


0.092±0.012 


0.13±0.02 






1.50 


2.00 


0.001±0.003 


0.018±0.004 


0.059±0.009 


0.13±0.02 






2.00 


2.50 




0.012±0.003 


0.036±0.006 


0.074±0.011 






2.50 


3.00 


± 


0.004±0.002 


0.025±0.005 


0.051±0.009 






3.00 


3.50 


± 


0.001±0.001 


0.008±0.002 


0.023±0.005 






3.50 


4.00 


± 


*_|_* 


0.005±0.002 


0.011±0.003 






4.00 


5.00 


± 


*_|_* 




0.007±0.002 






5.00 


6.50 


± 


± 




*_|_* 






6.50 


8.00 


± 


± 


± 


=f=_|_* 


0.200 


0.250 


0.50 


1.00 


0.05±0.03 


0.054±0.014 


0.13±0.02 


0.15±0.02 






1.00 


1.50 


0.004±0.007 


0.041±0.010 


0.096±0.014 


0.12±0.02 






1.50 


2.00 


*_L_* 


0.011±0.005 


0.08±0.02 


0.11±0.02 






2.00 


2.50 




0.007±0.003 


0.024±0.007 


0.060±0.014 






2.50 


3.00 


± 


0.001±0.001 


0.009±0.003 


0.024±0.007 






3.00 


3.50 


± 


0.001±0.001 


0.003±0.002 


0.010±0.004 






3.50 


4.00 


± 




0.001±0.001 


0.005±0.003 






4.00 


5.00 


± 




0.000±0.001 


O.OOliO.001 






5.00 


6.50 


± 


± 


0.000±0.001 


O.OOliO.001 






6.50 


8.00 


± 


± 


± 


O.OOOiO.001 
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TABLE XII: HARP results for the double-differential n production cross-section in the laboratory system, 
d 2 a n / (dpdQ.) , for p — N2 interactions at 12 GeV/c. Each row refers to a different (p mm <p< pmax, # m in < 6 < 
#max) bin, where p and 6 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 


d 2 a 7T+ /(dpdQ) 


d 2 a" /(dpdQ) 


(rad) 


(rad) (GeV/c) 


(GeV/c) 


(barn/(sr GeV/c)) (barn/(sr GeV/c)) 










12 GeV/c 


12 GeV/c 


0.05 


0.10 


0.50 


1.00 


0.24±0.04 


0.20±0.03 






1.00 


1.50 


0.32±0.03 


0.26±0.03 






1.50 


2.00 


0.29±0.03 


0.20±0.02 






2.00 


2.50 


0.31±0.03 


0.17±0.02 






2.50 


3.00 


0.28±0.02 


0.127±0.015 






3.00 


3.50 


0.18±0.02 


0.071±0.010 






3.50 


4.00 


0.15±0.02 


0.065±0.009 






4.00 


5.00 


0.098±0.009 


0.037±0.006 






5.00 


6.50 


0.046±0.006 


0.014±0.003 






6.50 


8.00 


0.014±0.003 


0.002±0.001 


0.10 


0.15 


0.50 


1.00 


0.24±0.04 


0.30±0.04 






1.00 


1.50 


0.34±0.03 


0.22±0.02 






1.50 


2.00 


0.25±0.03 


0.18±0.02 






2.00 


2.50 


0.23±0.03 


0.08±0.02 






2.50 


3.00 


0.15±0.02 


0.079±0.012 






3.00 


3.50 


0.102±0.013 


0.063±0.009 






3.50 


4.00 


0.062±0.010 


0.039±0.006 






4.00 


5.00 


0.034±0.005 


0.014±0.003 






5.00 


6.50 


0.011±0.002 


0.004±0.001 






6.50 


8.00 


0.002±0.001 




0.15 


0.20 


0.50 


1.00 


0.34±0.05 


0.27±0.04 






1.00 


1.50 


0.24±0.03 


0.15±0.02 






1.50 


2.00 


0.15±0.02 


0.12±0.02 






2.00 


2.50 


0.11±0.02 


0.084±0.015 






2.50 


3.00 


0.076±0.014 


0.042±0.011 






3.00 


3.50 


0.036±0.008 


0.017±0.007 






3.50 


4.00 


0.028±0.005 


0.020±0.005 






4.00 


5.00 


0.018±0.003 


0.006±0.002 






5.00 


6.50 


0.005±0.002 


± 






6.50 


8.00 


0.001±0.000 


± 


0.20 


0.25 


0.50 


1.00 


0.21±0.04 


0.15±0.03 






1.00 


1.50 


0.08±0.02 


0.19±0.03 






1.50 


2.00 


0.10±0.02 


0.12±0.03 






2.00 


2.50 


0.063±0.014 


0.08±0.02 






2.50 


3.00 


0.048±0.011 


0.050±0.012 






3.00 


3.50 


0.022±0.006 


0.015±0.005 






3.50 


4.00 


0.009±0.004 


0.006±0.004 






4.00 


5.00 


0.008±0.003 


0.001±0.002 






5.00 


6.50 


0.002±0.001 


O.OOOiO.001 






6.50 


8.00 


0.001±0.002 


O.OOOiO.001 
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TABLE XIII: HARP results for the double-differential n production cross-section in the laboratory system, 
d 2 cr 7r / (dpdQ.) , for p — O2 interactions at 12 GeV/c. Each row refers to a different (p m j n <p< pmax, # m in < 6 < 
#max) bin, where p and 6 are the pion momentum and polar angle, respectively. The central value as well as the 
square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


pmax 


d 2 a" + /(dpdQ.) 


cPa" /(dpdQ) 


(rad) 


(rad) (GeV/c) 


(GeV/c) 


(barn/(sr GeV/c) (barn/(sr GeV/c) 










12 GeV/c 


12 GeV/c 


0.05 


0.10 


0.50 


1.00 


0.29±0.06 


0.29±0.06 






1.00 


1.50 


0.42±0.06 


0.31±0.05 






1.50 


2.00 


0.44±0.05 


0.24±0.04 






2.00 


2.50 


0.37±0.05 


0.17±0.03 






2.50 


3.00 


0.31±0.04 


0.14±0.03 






3.00 


3.50 


0.27±0.04 


0.12±0.02 






3.50 


4.00 


0.18±0.03 


0.07±0.02 






4.00 


5.00 


0.13±0.02 


0.042±0.010 






5.00 


6.50 


0.063±0.010 


0.021±0.006 






6.50 


8.00 


0.025±0.007 


0.007±0.003 


0.10 


0.15 


0.50 


1.00 


0.33±0.06 


0.46±0.08 






1.00 


1.50 


0.49±0.07 


0.33±0.05 






1.50 


2.00 


0.32±0.05 


0.24±0.04 






2.00 


2.50 


0.33±0.05 


0.11±0.03 






2.50 


3.00 


0.22±0.04 


0.14±0.03 






3.00 


3.50 


0.11±0.02 


0.07±0.02 






3.50 


4.00 


0.10±0.02 


0.06±0.02 






4.00 


5.00 


0.052±0.011 


0.024±0.008 






5.00 


6.50 


0.015±0.005 


0.004±0.002 






6.50 


8.00 


0.002±0.001 


O.OOliO.001 


0.15 


0.20 


0.50 


1.00 


0.34±0.07 


0.24±0.06 






1.00 


1.50 


0.31±0.05 


0.20±0.04 






1.50 


2.00 


0.23±0.04 


0.11±0.03 






2.00 


2.50 


0.14±0.03 


0.10±0.03 






2.50 


3.00 


0.08±0.02 


0.05±0.02 






3.00 


3.50 


0.06±0.02 


0.015±0.010 






3.50 


4.00 


0.029±0.011 


0.027±0.013 






4.00 


5.00 


0.020±0.008 


0.011±0.007 






5.00 


6.50 


0.006±0.004 


0.001±0.002 






6.50 


8.00 


O.OOliO.001 


O.OOOiO.001 


0.20 


0.25 


0.50 


1.00 


0.23±0.06 


0.19±0.05 






1.00 


1.50 


0.14±0.04 


0.20±0.06 






1.50 


2.00 


0.13±0.04 


0.12±0.04 






2.00 


2.50 


0.07±0.03 


0.12±0.04 






2.50 


3.00 


0.04±0.02 


0.06±0.03 






3.00 


3.50 


0.04±0.02 


0.022±0.015 






3.50 


4.00 


0.018±0.013 


0.015±0.015 






4.00 


5.00 


0.011±0.009 


0.005±0.006 






5.00 


6.50 


0.003±0.005 


0.002±0.007 






6.50 


8.00 


0.002±0.005 


0.001±0.007 



38 



TABLE XIV: HARP results for the double-differential ix + production cross-section in the laboratory system, 
d 2 er 7r / (dpdQ) , for p-Al interactions at 3, 5, 8, 12, 12.9 GeV/c. Each row refers to a different (p mm < p < 
Pma.x,8 m [ n < < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central 
value as well as the square-root of the diagonal elements of the covariance matrix are given. 



min 


#max 


Pmin 


JJmax 






d 2 a n+ /(dpdQ) 




(rad) 


(rad) (GeV/c) 


(GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.19±0.06 


0.17±0.03 


0.34±0.04 


0.47±0.06 


0.49±0.04 






1.00 


1.50 


0.02±0.02 


0.27±0.03 


0.35±0.03 


0.42±0.04 


0.46±0.03 






1.50 


2.00 




0.20±0.03 


0.36±0.03 


0.48±0.04 


0.52±0.03 






2.00 


2.50 




0.064±0.012 


0.31±0.02 


0.47±0.04 


0.47±0.02 






2.50 


3.00 


± 


0.038±0.010 


0.23±0.02 


0.44±0.04 


0.41±0.02 






3.00 


3.50 


± 


0.006±0.005 


0.18±0.02 


0.28±0.03 


0.36±0.02 






3.50 


4.00 


± 


0.004±0.003 


0.109±0.014 


0.23±0.03 


0.27±0.02 






4.00 


5.00 


± 


0.001±0.001 


0.077±0.010 


0.143±0.014 


0.169±0.011 






5.00 


6.50 


± 


± 


0.006±0.002 


0.048±0.007 


0.077±0.005 






6.50 


8.00 


± 


± 


± 


0.021±0.004 


0.029±0.003 


0.100 


0.150 


0.50 


1.00 


0.13±0.05 


0.24±0.04 


0.36±0.04 


0.44±0.06 


0.49±0.04 






1.00 


1.50 


0.01±0.01 


0.22±0.03 


0.38±0.04 


0.50±0.05 


0.53±0.04 






1.50 


2.00 




0.15±0.02 


0.32±0.03 


0.46±0.05 


0.47±0.03 






2.00 


2.50 




0.048±0.011 


0.19±0.02 


0.32±0.04 


0.35±0.03 






2.50 


3.00 


± 


0.015±0.007 


0.13±0.02 


0.22±0.03 


0.28±0.02 






3.00 


3.50 


± 


0.003±0.001 


0.091±0.013 


0.18±0.02 


0.18±0.02 






3.50 


4.00 


± 


0.001±0.001 


0.053±0.008 


0.09±0.02 


0.128±0.011 






4.00 


5.00 


± 




0.020±0.005 


0.050±0.008 


0.070±0.007 






5.00 


6.50 


± 


± 


0.001±0.000 


0.017±0.003 


0.021±0.003 






6.50 


8.00 


± 


± 


± 


0.004±0.001 


0.004±0.001 


0.150 


0.200 


0.50 


1.00 


0.19±0.07 


0.26±0.04 


0.39±0.05 


0.40±0.06 


0.50±0.05 






1.00 


1.50 


0.03±0.02 


0.17±0.03 


0.31±0.03 


0.34±0.04 


0.35±0.03 






1.50 


2.00 


*_L.* 


0.054±0.013 


0.17±0.02 


0.29±0.03 


0.26±0.02 






2.00 


2.50 


*_|_* 


0.023±0.010 


0.098±0.014 


0.17±0.02 


0.17±0.02 






2.50 


3.00 


± 


0.005±0.004 


0.058±0.010 


0.10±0.02 


0.117±0.011 






3.00 


3.50 


± 


0.001±0.001 


0.039±0.007 


0.046±0.011 


0.072±0.008 






3.50 


4.00 


± 


*_|_* 


0.023±0.006 


0.033±0.008 


0.047±0.006 






4.00 


5.00 


± 


*_|_* 


0.008±0.002 


0.025±0.007 


0.021±0.004 






5.00 


6.50 


± 


± 




0.005±0.003 


0.005±0.001 






6.50 


8.00 


± 


± 


± 


*_|_=f= 




0.200 


0.250 


0.50 


1.00 


0.05±0.04 


0.16±0.03 


0.24±0.03 


0.27±0.05 


0.33±0.04 






1.00 


1.50 


0.01±0.02 


0.08±0.02 


0.18±0.03 


0.17±0.04 


0.21±0.02 






1.50 


2.00 


*_|_=f= 


0.039±0.014 


0.17±0.03 


0.10±0.03 


0.15±0.02 






2.00 


2.50 


*_|_* 


0.020±0.009 


0.07±0.02 


0.10±0.03 


0.103±0.014 






2.50 


3.00 


± 


0.004±0.003 


0.037±0.010 


0.05±0.02 


0.076±0.010 






3.00 


3.50 


± 


0.001±0.001 


0.025±0.007 


0.031±0.010 


0.038±0.006 






3.50 


4.00 


± 


0.001±0.001 


0.013±0.005 


0.018±0.009 


0.021±0.004 






4.00 


5.00 


± 


0.001±0.002 


0.007±0.004 


0.012±0.008 


0.015±0.003 






5.00 


6.50 


± 


± 


0.001±0.002 


0.005±0.004 


0.004±0.001 






6.50 


8.00 


± 


± 


± 


0.002±0.007 


0.001±0.001 
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TABLE XV: HARP results for the double-differential n~ production cross-section in the laboratory system, 
d 2 er 7r / (dpdQ,) , for p-Al interactions at 3, 5, 8, 12, 12.9 GeV/c. Each row refers to a different (p mm < p < 
Pma.x,8 m [ n < < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central 
value as well as the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


JJmax 




d 2 a n /(dpdQ,) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.03±0.03 


0.14±0.03 


0.28±0.04 


0.38±0.05 


0.37±0.04 






1.00 


1.50 




0.10±0.02 


0.21±0.02 


0.35±0.04 


0.40±0.03 






1.50 


2.00 


0.001±0.003 


0.050±0.013 


0.19±0.02 


0.34±0.04 


0.36±0.02 






2.00 


2.50 




0.023±0.008 


0.117±0.014 


0.27±0.03 


0.29±0.02 






2.50 


3.00 


± 


0.007±0.004 


0.066±0.013 


0.17±0.02 


0.20±0.02 






3.00 


3.50 


± 




0.040±0.007 


0.14±0.02 


0.179±0.011 






3.50 


4.00 


± 




0.032±0.006 


0.11±0.02 


0.115±0.010 






4.00 


5.00 


± 


±* 


0.011±0.004 


0.048±0.009 


0.069±0.006 






5.00 


6.50 


± 


± 


O.OOliO.001 


0.020±0.004 


0.028±0.003 






6.50 


8.00 


± 


± 


± 


0.003±0.002 


0.007±0.002 


0.100 


0.150 


0.50 


1.00 


0.03±0.03 


0.16±0.03 


0.32±0.04 


0.48±0.07 


0.46±0.05 






1.00 


1.50 




0.08±0.02 


0.22±0.02 


0.39±0.04 


0.38±0.03 






1.50 


2.00 


0.002±0.005 


0.038±0.010 


0.14±0.02 


0.24±0.03 


0.30±0.02 






2.00 


2.50 




0.017±0.006 


0.080±0.012 


0.19±0.03 


0.23±0.02 






2.50 


3.00 


± 


0.011±0.004 


0.038±0.007 


0.14±0.02 


0.172±0.015 






3.00 


3.50 


± 


0.002±0.001 


0.026±0.006 


0.095±0.015 


0.103±0.009 






3.50 


4.00 


± 




0.014±0.003 


0.066±0.012 


0.081±0.008 






4.00 


5.00 


± 




0.003±0.002 


0.032±0.007 


0.038±0.005 






5.00 


6.50 


± 


± 




0.008±0.002 


0.007±0.002 






6.50 


8.00 


± 


± 


± 


0.001±0.000 


0.001±0.000 


0.150 


0.200 


0.50 


1.00 


0.07±0.04 


0.18±0.03 


0.35±0.04 


0.49±0.07 


0.43±0.04 






1.00 


1.50 


0.02±0.02 


0.08±0.02 


0.17±0.02 


0.25±0.04 


0.35±0.03 






1.50 


2.00 


*_|_* 


0.027±0.008 


0.093±0.014 


0.19±0.03 


0.22±0.02 






2.00 


2.50 


*_|_* 


0.014±0.005 


0.065±0.011 


0.12±0.02 


0.156±0.015 






2.50 


3.00 


± 


0.005±0.004 


0.040±0.008 


0.064±0.013 


0.084±0.010 






3.00 


3.50 


± 


0.000±0.001 


0.013±0.004 


0.063±0.014 


0.046±0.006 






3.50 


4.00 


± 




0.005±0.002 


0.022±0.007 


0.026±0.004 






4.00 


5.00 


± 




O.OOliO.001 


0.007±0.003 


0.012±0.002 






5.00 


6.50 


± 


± 




0.001±0.001 


O.OOliO.001 






6.50 


8.00 


± 


± 


± 


*_|_=f= 




0.200 


0.250 


0.50 


1.00 


0.03±0.03 


0.12±0.03 


0.21±0.03 


0.29±0.05 


0.34±0.04 






1.00 


1.50 


0.02±0.03 


0.07±0.02 


0.19±0.03 


0.26±0.04 


0.28±0.03 






1.50 


2.00 


O.OOOiO.001 


0.010±0.007 


0.14±0.03 


0.16±0.04 


0.26±0.04 






2.00 


2.50 




0.005±0.004 


0.06±0.02 


0.11±0.03 


0.13±0.02 






2.50 


3.00 


± 


0.002±0.002 


0.022±0.007 


0.028±0.012 


0.064±0.012 






3.00 


3.50 


± 




0.009±0.004 


0.016±0.007 


0.031±0.006 






3.50 


4.00 


± 




0.002±0.002 


0.013±0.006 


0.013±0.003 






4.00 


5.00 


± 




0.001±0.003 


0.004±0.004 


0.005±0.001 






5.00 


6.50 


± 


± 


0.000±0.001 


0.000±0.001 


O.OOliO.001 






6.50 


8.00 


± 


± 


± 


0.000±0.004 


O.OOliO.001 
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TABLE XVI: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 er 7r / (dpdfi) , for p-Cu interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < p < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 




d 2 a' K+ /(dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.33±0.13 


0.30±0.06 


0.55±0.07 


0.77±0.10 






1.00 


1.50 


0.05±0.04 


0.40±0.05 


0.59±0.05 


0.76±0.07 






1.50 


2.00 


0.02±0.02 


0.25±0.04 


0.56±0.04 


0.81±0.07 






2.00 


2.50 




0.10±0.02 


0.44±0.04 


0.75±0.06 






2.50 


3.00 


± 


0.025±0.011 


0.26±0.03 


0.61±0.05 






3.00 


3.50 


± 


0.025±0.013 


0.13±0.02 


0.53±0.05 






3.50 


4.00 


± 


0.021±0.012 


0.05±0.02 


0.36±0.04 






4.00 


5.00 


± 


0.001±0.002 


0.036±0.010 


0.22±0.02 






5.00 


6.50 


± 


± 


0.005±0.002 


0.087±0.012 






6.50 


8.00 


± 


± 


± 


0.028±0.006 


0.100 


0.150 


0.50 


1.00 


0.19±0.10 


0.49±0.07 


0.63±0.07 


0.90±0.10 






1.00 


1.50 


0.07±0.06 


0.36±0.05 


0.60±0.06 


0.84±0.08 






1.50 


2.00 


0.002±0.006 


0.20±0.03 


0.48±0.04 


0.72±0.07 






2.00 


2.50 




0.05±0.02 


0.28±0.03 


0.53±0.06 






2.50 


3.00 


± 


0.012±0.007 


0.12±0.02 


0.37±0.04 






3.00 


3.50 


± 


0.007±0.003 


0.057±0.013 


0.27±0.03 






3.50 


4.00 


± 


0.001±0.001 


0.022±0.006 


0.18±0.02 






4.00 


5.00 


± 




0.014±0.004 


0.09±0.02 






5.00 


6.50 


± 


± 


0.002±0.001 


0.015±0.004 






6.50 


8.00 


± 


± 


± 


0.003±0.001 


0.150 


0.200 


0.50 


1.00 


0.26±0.13 


0.47±0.07 


0.61±0.07 


0.86±0.11 






1.00 


1.50 


0.05±0.05 


0.25±0.04 


0.42±0.04 


0.64±0.07 






1.50 


2.00 


0.00±0.01 


0.09±0.02 


0.28±0.03 


0.39±0.05 






2.00 


2.50 




0.016±0.011 


0.15±0.02 


0.31±0.04 






2.50 


3.00 


± 


0.006±0.006 


0.076±0.015 


0.14±0.03 






3.00 


3.50 


± 


0.004±0.003 


0.040±0.010 


0.13±0.02 






3.50 


4.00 


± 


0.000±0.001 


0.009±0.004 


0.07±0.02 






4.00 


5.00 


± 




0.010±0.004 


0.019±0.006 






5.00 


6.50 


± 


± 


0.001±0.001 


0.005±0.002 






6.50 


8.00 


± 


± 


± 


0.002±0.001 


0.200 


0.250 


0.50 


1.00 


0.04±0.07 


0.32±0.06 


0.47±0.06 


0.59±0.09 






1.00 


1.50 


0.03±0.13 


0.14±0.04 


0.24±0.04 


0.33±0.06 






1.50 


2.00 




0.04±0.02 


0.21±0.04 


0.14±0.03 






2.00 


2.50 




0.04±0.02 


0.09±0.02 


0.12±0.03 






2.50 


3.00 


± 


0.002±0.004 


0.028±0.011 


0.11±0.03 






3.00 


3.50 


± 


0.001±0.001 


0.018±0.008 


0.03±0.02 






3.50 


4.00 


± 


0.001±0.002 


0.017±0.007 


O.OlOiO.011 






4.00 


5.00 


± 


0.000±0.001 


0.015±0.007 


0.009±0.010 






5.00 


6.50 


± 


± 


0.002±0.003 


0.001±0.002 






6.50 


8.00 


± 


± 


± 


0.001±0.009 
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TABLE XVII: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 cr 7r / (dpdQ) , for p-Cu interactions at 3, 5, 8, 12 GcV/c. Each row refers to a different (p mm < p < Pmax, "min — 
9 < 9max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 




d 2 a n /(dpdQ,) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.07±0.06 


0.18±0.04 


0.45±0.06 


0.58±0.08 






1.00 


1.50 




0.18±0.03 


0.38±0.04 


0.70±0.07 






1.50 


2.00 


0.002±0.008 


0.08±0.02 


0.32±0.03 


0.57±0.06 






2.00 


2.50 




0.05±0.02 


0.17±0.02 


0.41±0.05 






2.50 


3.00 


± 


0.003±0.003 


0.08±0.02 


0.35±0.04 






3.00 


3.50 


± 


0.001±0.002 


0.041±0.009 


0.22±0.04 






3.50 


4.00 


± 


*_|_* 


0.050±0.010 


0.15±0.02 






4.00 


5.00 


± 


=f=_|_* 


0.012±0.006 


0.089±0.014 






5.00 


6.50 


± 


± 


0.000±0.001 


0.029±0.007 






6.50 


8.00 


± 


± 


± 


0.006±0.004 


0.100 


0.150 


0.50 


1.00 


0.14±0.10 


0.26±0.05 


0.57±0.07 


0.90±0.12 






1.00 


1.50 




0.13±0.03 


0.42±0.04 


0.71±0.07 






1.50 


2.00 




0.053±0.015 


0.22±0.03 


0.38±0.05 






2.00 


2.50 




0.020±0.009 


0.14±0.02 


0.31±0.04 






2.50 


3.00 


± 


0.008±0.004 


0.088±0.014 


0.22±0.03 






3.00 


3.50 


± 


0.001±0.001 


0.044±0.009 


0.14±0.02 






3.50 


4.00 


± 


*_|_* 


0.020±0.006 


0.10±0.02 






4.00 


5.00 


± 


*_|_* 


0.006±0.002 


0.033±0.008 






5.00 


6.50 


± 


± 




0.007±0.003 






6.50 


8.00 


± 


± 


± 


0.001±0.000 


0.150 


0.200 


0.50 


1.00 


0.04±0.04 


0.26±0.05 


0.60±0.07 


0.95±0.12 






1.00 


1.50 


0.02±0.03 


0.11±0.03 


0.31±0.04 


0.42±0.06 






1.50 


2.00 


0.01±0.03 


0.05±0.02 


0.19±0.03 


0.40±0.05 






2.00 


2.50 


0.000±0.002 


0.027±0.010 


0.10±0.02 


0.22±0.04 






2.50 


3.00 


± 


0.018±0.010 


0.048±0.011 


0.13±0.03 






3.00 


3.50 


± 


0.001±0.001 


0.021±0.007 


0.08±0.02 






3.50 


4.00 


± 


*_|_* 


0.010±0.004 


0.048±0.014 






4.00 


5.00 


± 


*_|_* 


0.002±0.001 


0.029±0.008 






5.00 


6.50 


± 


± 




0.003±0.002 






6.50 


8.00 


± 


± 


± 


*_|_=f= 


0.200 


0.250 


0.50 


1.00 


0.09±0.10 


0.12±0.04 


0.42±0.06 


0.56±0.09 






1.00 


1.50 




0.11±0.03 


0.27±0.04 


0.34±0.06 






1.50 


2.00 




0.03±0.02 


0.19±0.04 


0.37±0.07 






2.00 


2.50 




0.006±0.005 


0.051±0.015 


0.23±0.05 






2.50 


3.00 


± 


0.001±0.001 


0.016±0.007 


0.08±0.02 






3.00 


3.50 


± 


0.000±0.001 


0.006±0.005 


0.033±0.012 






3.50 


4.00 


± 


0.000±0.002 


0.003±0.003 


0.026±0.013 






4.00 


5.00 


± 


0.000±0.002 


0.001±0.003 


0.011±0.009 






5.00 


6.50 


± 


± 


0.000±0.001 


0.005±0.008 






6.50 


8.00 


± 


± 


± 


0.001±0.005 
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TABLE XVIII: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 er 7r / (dpdfl) , for p-Sn interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (f> mm < p < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 




d 2 a' K+ /(dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.32±0.12 


0.42±0.08 


0.80±0.10 


1.18±0.12 






1.00 


1.50 


0.001±0.003 


0.39±0.06 


0.83±0.07 


1.07±0.09 






1.50 


2.00 


0.016±0.011 


0.25±0.05 


0.74±0.06 


1.09±0.07 






2.00 


2.50 


0.01±0.01 


0.10±0.03 


0.68±0.06 


1.02±0.07 






2.50 


3.00 


± 


0.05±0.02 


0.38±0.05 


0.88±0.06 






3.00 


3.50 


± 


*_|_* 


0.14±0.03 


0.71±0.05 






3.50 


4.00 


± 


*_|_* 


0.06±0.02 


0.44±0.05 






4.00 


5.00 


± 


=f=_|_* 


0.06±0.02 


0.31±0.02 






5.00 


6.50 


± 


± 


0.008±0.003 


0.12±0.02 






6.50 


8.00 


± 


± 


± 


0.044±0.007 


0.100 


0.150 


0.50 


1.00 


0.10±0.07 


0.53±0.08 


1.03±0.11 


1.24±0.12 






1.00 


1.50 


0.04±0.05 


0.45±0.07 


0.71±0.07 


1.21±0.10 






1.50 


2.00 


0.02±0.02 


0.21±0.04 


0.66±0.06 


0.94±0.08 






2.00 


2.50 




0.09±0.03 


0.39±0.04 


0.75±0.07 






2.50 


3.00 


± 


0.05±0.02 


0.18±0.03 


0.60±0.05 






3.00 


3.50 


± 


0.008±0.004 


0.054±0.015 


0.35±0.04 






3.50 


4.00 


± 


0.002±0.001 


0.031±0.009 


0.20±0.02 






4.00 


5.00 


± 


0.001±0.001 


0.016±0.005 


0.116±0.015 






5.00 


6.50 


± 


± 


0.002±0.001 


0.029±0.006 






6.50 


8.00 


± 


± 


± 


0.005±0.001 


0.150 


0.200 


0.50 


1.00 


0.19±0.12 


0.43±0.08 


0.91±0.10 


1.16±0.13 






1.00 


1.50 


0.03±0.05 


0.30±0.05 


0.55±0.06 


1.01±0.09 






1.50 


2.00 


=t=_|_* 


0.11±0.03 


0.38±0.04 


0.60±0.06 






2.00 


2.50 


0.01±0.05 


0.04±0.02 


0.20±0.03 


0.38±0.04 






2.50 


3.00 


± 


0.01±0.01 


0.09±0.02 


0.22±0.03 






3.00 


3.50 


± 


0.000±0.001 


0.043±0.011 


0.15±0.02 






3.50 


4.00 


± 


0.002±0.002 


0.019±0.007 


0.089±0.014 






4.00 


5.00 


± 


0.000±0.001 


0.007±0.003 


0.042±0.009 






5.00 


6.50 


± 


± 


0.001±0.001 


0.014±0.004 






6.50 


8.00 


± 


± 


± 


O.OOliO.001 


0.200 


0.250 


0.50 


1.00 


*_j_* 


0.39±0.08 


0.73±0.10 


0.85±0.11 






1.00 


1.50 


*_|_* 


0.16±0.05 


0.48±0.07 


0.50±0.07 






1.50 


2.00 




0.07±0.03 


0.30±0.05 


0.32±0.05 






2.00 


2.50 




0.015±0.010 


0.22±0.04 


0.21±0.04 






2.50 


3.00 


± 


0.012±0.010 


0.08±0.02 


0.16±0.03 






3.00 


3.50 


± 


0.013±0.009 


0.033±0.011 


0.08±0.02 






3.50 


4.00 


± 


0.006±0.006 


0.018±0.008 


0.039±0.011 






4.00 


5.00 


± 


0.002±0.003 


0.013±0.007 


0.028±0.009 






5.00 


6.50 


± 


± 


0.004±0.004 


0.006±0.003 






6.50 


8.00 


± 


± 


± 


0.003±0.006 
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TABLE XIX: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 a 7 ' I (dpdfl) , for p-Sn interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (f> mm < V < Pmax, "min — 
9 < 9max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


^min 


Pmax 




d 2 a n /{dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.11±0.07 


0.26±0.07 


0.60±0.08 


0.92±0.10 






1.00 


1.50 


0.004±0.006 


0.21±0.04 


0.50±0.05 


1.02±0.08 






1.50 


2.00 


0.01±0.01 


0.15±0.04 


0.41±0.04 


0.82±0.06 






2.00 


2.50 




0.033±0.015 


0.24±0.03 


0.50±0.05 






2.50 


3.00 


± 


0.012±0.008 


0.13±0.02 


0.46±0.04 






3.00 


3.50 


± 


*_|_* 


0.09±0.02 


0.30±0.03 






3.50 


4.00 


± 


*_|_* 


0.044±0.010 


0.20±0.02 






4.00 


5.00 


± 


*_|_* 


0.024±0.009 


0.12±0.02 






5.00 


6.50 


± 


± 


0.003±0.002 


0.030±0.007 






6.50 


8.00 


± 


± 


± 


0.006±0.003 


0.100 


0.150 


0.50 


1.00 


0.10±0.07 


0.36±0.08 


0.89±0.11 


1.40±0.15 






1.00 


1.50 




0.11±0.03 


0.59±0.06 


0.92±0.08 






1.50 


2.00 




0.07±0.02 


0.32±0.04 


0.63±0.06 






2.00 


2.50 




0.024±0.012 


0.18±0.03 


0.43±0.04 






2.50 


3.00 


± 


0.016±0.007 


0.11±0.02 


0.25±0.03 






3.00 


3.50 


± 


0.007±0.005 


0.064±0.014 


0.22±0.03 






3.50 


4.00 


± 


*_|_* 


0.020±0.006 


0.12±0.02 






4.00 


5.00 


± 


*_|_* 


0.010±0.004 


0.058±0.011 






5.00 


6.50 


± 


± 


0.001±0.001 


0.008±0.003 






6.50 


8.00 


± 


± 


± 


0.001±0.000 


0.150 


0.200 


0.50 


1.00 


0.06±0.11 


0.34±0.07 


0.77±0.09 


1.17±0.13 






1.00 


1.50 


0.06±0.04 


0.15±0.04 


0.49±0.06 


0.75±0.07 






1.50 


2.00 




0.07±0.02 


0.21±0.03 


0.51±0.05 






2.00 


2.50 


0.005±0.012 


0.06±0.02 


0.14±0.02 


0.30±0.04 






2.50 


3.00 


± 


0.00±0.01 


0.07±0.02 


0.18±0.03 






3.00 


3.50 


± 


=f=_|_* 


0.016±0.008 


0.10±0.02 






3.50 


4.00 


± 


*_|_* 


0.012±0.005 


0.061±0.013 






4.00 


5.00 


± 


*_|_* 


0.001±0.001 


0.033±0.007 






5.00 


6.50 


± 


± 




0.005±0.002 






6.50 


8.00 


± 


± 


± 




0.200 


0.250 


0.50 


1.00 


0.05±0.06 


0.36±0.08 


0.59±0.08 


0.88±0.11 






1.00 


1.50 


0.05±0.05 


0.12±0.04 


0.42±0.06 


0.54±0.07 






1.50 


2.00 


*_L_* 


0.02±0.02 


0.29±0.05 


0.55±0.09 






2.00 


2.50 




0.005±0.005 


0.16±0.03 


0.27±0.05 






2.50 


3.00 


± 


*_|_=f= 


0.06±0.02 


0.13±0.03 






3.00 


3.50 


± 


0.000±0.001 


0.023±0.011 


0.06±0.02 






3.50 


4.00 


± 




0.004±0.004 


0.033±0.011 






4.00 


5.00 


± 


*_l_=f= 


0.001±0.005 


0.008±0.004 






5.00 


6.50 


± 


± 


0.000±0.001 


0.002±0.003 






6.50 


8.00 


± 


± 


± 


0.000±0.003 
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TABLE XX: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d 2 er 7r / (dpdfl) , for p-Ta interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < P < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


^min 


Pmax 




d 2 a w+ /(dpduj) 




(rad) 


(rad) (GeV/c) 


(GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.26±0.11 


0.55±0.12 


0.87±0.12 


1.65±0.20 






1.00 


1.50 


0.004±0.007 


0.54±0.10 


0.95±0.09 


1.27±0.13 






1.50 


2.00 


0.014±0.009 


0.25±0.07 


0.80±0.08 


1.37±0.12 






2.00 


2.50 


0.01±0.01 


0.11±0.03 


0.59±0.07 


1.05±0.11 






2.50 


3.00 


± 


0.05±0.03 


0.34±0.06 


1.19±0.11 






3.00 


3.50 


± 


*_|_* 


0.15±0.05 


0.67±0.08 






3.50 


4.00 


± 


*_|_* 


0.05±0.03 


0.56±0.07 






4.00 


5.00 


± 


*_|_* 


0.05±0.02 


0.30±0.04 






5.00 


6.50 


± 


± 


0.007±0.003 


0.12±0.02 






6.50 


8.00 


± 


± 


± 


0.038±0.009 


0.100 


0.150 


0.50 


1.00 


0.22±0.12 


0.56±0.11 


0.93±0.11 


1.60±0.19 






1.00 


1.50 


0.03±0.04 


0.39±0.09 


0.82±0.10 


1.60±0.16 






1.50 


2.00 


0.002±0.008 


0.26±0.07 


0.82±0.09 


1.22±0.13 






2.00 


2.50 


0.002±0.007 


0.04±0.02 


0.40±0.06 


0.84±0.10 






2.50 


3.00 


± 


0.02±0.02 


0.15±0.04 


0.61±0.07 






3.00 


3.50 


± 


0.009±0.009 


0.09±0.03 


0.41±0.05 






3.50 


4.00 


± 




0.02±0.02 


0.26±0.05 






4.00 


5.00 


± 


0.000±0.001 


0.009±0.004 


0.15±0.02 






5.00 


6.50 


± 


± 


0.002±0.001 


0.034±0.009 






6.50 


8.00 


± 


± 


± 


0.009±0.003 


0.150 


0.200 


0.50 


1.00 


0.38±0.17 


0.57±0.12 


0.93±0.12 


1.37±0.19 






1.00 


1.50 


0.13±0.09 


0.39±0.08 


0.77±0.09 


1.15±0.13 






1.50 


2.00 




0.10±0.04 


0.38±0.06 


0.66±0.09 






2.00 


2.50 




0.06±0.03 


0.20±0.04 


0.48±0.08 






2.50 


3.00 


± 


O.OOiO.OO 


0.12±0.03 


0.34±0.06 






3.00 


3.50 


± 


O.OOiO.OO 


0.06±0.02 


0.12±0.03 






3.50 


4.00 


± 


0.000±0.001 


0.016±0.009 


0.08±0.02 






4.00 


5.00 


± 




0.007±0.004 


0.06±0.02 






5.00 


6.50 


± 


± 


0.002±0.002 


0.011±0.005 






6.50 


8.00 


± 


± 


± 


0.003±0.002 


0.200 


0.250 


0.50 


1.00 


0.16±0.34 


0.48±0.11 


0.63±0.10 


1.24±0.18 






1.00 


1.50 


0.04±0.38 


0.23±0.07 


0.46±0.08 


0.65±0.11 






1.50 


2.00 


*_L_* 


0.04±0.02 


0.20±0.05 


0.60±0.11 






2.00 


2.50 




0.06±0.03 


0.24±0.05 


0.49±0.10 






2.50 


3.00 


± 


0.03±0.02 


0.09±0.03 


0.24±0.06 






3.00 


3.50 


± 


0.020±0.015 


0.027±0.015 


0.10±0.03 






3.50 


4.00 


± 


0.007±0.008 


0.016±0.013 


0.04±0.02 






4.00 


5.00 


± 


0.005±0.008 


0.011±0.010 


0.05±0.03 






5.00 


6.50 


± 


± 


0.002±0.005 


0.015±0.012 






6.50 


8.00 


± 


± 


± 


0.003±0.017 
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TABLE XXI: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 a 7 ' I (dpdfl) , for p-Ta interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < P < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


^min 


Pmax 




d 2 a n /{dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.11±0.07 


0.31±0.10 


0.68±0.10 


1.21±0.16 






1.00 


1.50 


0.004±0.007 


0.29±0.07 


0.60±0.08 


1.23±0.12 






1.50 


2.00 


*_j_* 


0.09±0.04 


0.42±0.06 


0.89±0.10 






2.00 


2.50 


0.03±0.03 


0.06±0.03 


0.27±0.04 


0.92±0.10 






2.50 


3.00 


± 


0.015±0.012 


0.12±0.03 


0.47±0.08 






3.00 


3.50 


± 


0.007±0.007 


0.046±0.015 


0.42±0.06 






3.50 


4.00 


± 




0.09±0.02 


0.24±0.05 






4.00 


5.00 


± 


0.000±0.001 


0.017±0.011 


0.13±0.02 






5.00 


6.50 


± 


± 




0.044±0.012 






6.50 


8.00 


± 


± 


± 


0.010±0.006 


0.100 


0.150 


0.50 


1.00 


0.24±0.12 


0.30±0.09 


1.18±0.15 


1.79±0.22 






1.00 


1.50 




0.16±0.05 


0.68±0.08 


1.22±0.13 






1.50 


2.00 




0.08±0.03 


0.41±0.05 


0.71±0.09 






2.00 


2.50 




0.04±0.02 


0.17±0.03 


0.48±0.07 






2.50 


3.00 


± 


0.026±0.014 


0.10±0.02 


0.42±0.06 






3.00 


3.50 


± 


*_|_=f= 


0.07±0.02 


0.30±0.05 






3.50 


4.00 


± 


*_|_* 


0.031±0.010 


0.18±0.03 






4.00 


5.00 


± 


*_|_* 


0.014±0.006 


0.059±0.015 






5.00 


6.50 


± 


± 


0.001±0.001 


0.014±0.005 






6.50 


8.00 


± 


± 


± 


0.000±0.001 


0.150 


0.200 


0.50 


1.00 


0.14±0.09 


0.60±0.12 


0.86±0.12 


1.70±0.22 






1.00 


1.50 


0.05±0.05 


0.15±0.05 


0.55±0.07 


0.87±0.11 






1.50 


2.00 




0.07±0.02 


0.20±0.04 


0.68±0.09 






2.00 


2.50 




0.06±0.02 


0.16±0.03 


0.42±0.07 






2.50 


3.00 


± 


0.00±0.01 


0.06±0.02 


0.17±0.04 






3.00 


3.50 


± 


0.007±0.008 


0.031±0.014 


0.15±0.03 






3.50 


4.00 


± 






0.07±0.02 






4.00 


5.00 


± 






0.045±0.014 






5.00 


6.50 


± 


± 


0.000±0.001 


0.004±0.003 






6.50 


8.00 


± 


± 


± 


0.000±0.001 


0.200 


0.250 


0.50 


1.00 


0.14±0.11 


0.26±0.09 


0.76±0.11 


1.05±0.16 






1.00 


1.50 


0.04±0.05 


0.21±0.06 


0.32±0.06 


0.97±0.15 






1.50 


2.00 




0.03±0.03 


0.14±0.04 


0.75±0.14 






2.00 


2.50 




0.04±0.03 


0.06±0.02 


0.33±0.08 






2.50 


3.00 


± 


0.003±0.006 


0.021±0.010 


0.20±0.06 






3.00 


3.50 


± 


*_|_* 


*_L_* 


0.07±0.03 






3.50 


4.00 


± 


*_|_* 


0.001±0.002 


0.024±0.015 






4.00 


5.00 


± 


*_|_=f= 


0.000±0.005 


0.008±0.011 






5.00 


6.50 


± 


± 


0.000±0.002 


0.001±0.002 






6.50 


8.00 


± 


± 


± 


0.000±0.009 
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TABLE XXII: HARP results for the double-differential 7r + production cross-section in the laboratory system, 
d 2 a 7 ' I (dpdQ,) , for p-Pb interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < p < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 




d 2 a w+ /(dpduj) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.09±0.05 


0.45±0.11 


0.99±0.14 


1.78±0.27 






1.00 


1.50 




0.63±0.11 


0.86±0.09 


1.21±0.16 






1.50 


2.00 


0.012±0.009 


0.46±0.09 


0.81±0.08 


1.47±0.19 






2.00 


2.50 


0.01±0.01 


0.06±0.02 


0.66±0.08 


1.40±0.16 






2.50 


3.00 


± 


0.06±0.04 


0.41±0.07 


1.02±0.13 






3.00 


3.50 


± 


0.00±0.01 


0.16±0.05 


0.83±0.11 






3.50 


4.00 


± 


0.01±0.02 


0.03±0.06 


0.56±0.09 






4.00 


5.00 


± 


*_|_* 


0.06±0.03 


0.39±0.06 






5.00 


6.50 


± 


± 


0.01±0.01 


0.16±0.03 






6.50 


8.00 


± 


± 


± 


0.034±0.011 


0.100 


0.150 


0.50 


1.00 


0.35±0.17 


0.74±0.14 


1.18±0.14 


1.77±0.25 






1.00 


1.50 


0.04±0.05 


0.42±0.10 


0.77±0.10 


1.45±0.19 






1.50 


2.00 


0.01±0.02 


0.18±0.05 


0.80±0.08 


1.42±0.19 






2.00 


2.50 


0.003±0.007 


0.07±0.03 


0.45±0.07 


0.84±0.13 






2.50 


3.00 


± 


0.06±0.03 


0.15±0.04 


0.61±0.10 






3.00 


3.50 


± 


0.014±0.009 


0.05±0.02 


0.45±0.08 






3.50 


4.00 


± 


0.001±0.001 


0.031±0.012 


0.26±0.05 






4.00 


5.00 


± 




0.011±0.006 


0.13±0.03 






5.00 


6.50 


± 


± 


0.000±0.000 


0.044±0.013 






6.50 


8.00 


± 


± 


± 


0.008±0.004 


0.150 


0.200 


0.50 


1.00 


0.14±0.13 


0.66±0.14 


1.07±0.14 


1.32±0.23 






1.00 


1.50 


0.08±0.07 


0.30±0.08 


0.68±0.08 


1.08±0.16 






1.50 


2.00 


0.02±0.04 


0.14±0.05 


0.40±0.06 


0.75±0.12 






2.00 


2.50 


0.00±0.01 


0.04±0.04 


0.24±0.05 


0.56±0.10 






2.50 


3.00 


± 


0.01±0.01 


0.08±0.02 


0.39±0.08 






3.00 


3.50 


± 


0.000±0.001 


0.040±0.014 


0.12±0.04 






3.50 


4.00 


± 




0.018±0.008 


0.13±0.04 






4.00 


5.00 


± 




0.009±0.005 


0.05±0.02 






5.00 


6.50 


± 


± 


0.002±0.002 


0.015±0.008 






6.50 


8.00 


± 


± 


± 


0.015±0.010 


0.200 


0.250 


0.50 


1.00 


*_j_* 


0.49±0.13 


0.54±0.09 


0.97±0.20 






1.00 


1.50 


0.08±0.10 


0.21±0.08 


0.41±0.08 


0.61±0.15 






1.50 


2.00 


0.01±0.06 


0.06±0.04 


0.22±0.06 


0.46±0.13 






2.00 


2.50 




0.02±0.02 


0.13±0.03 


0.33±0.09 






2.50 


3.00 


± 


0.02±0.02 


0.05±0.02 


0.33±0.11 






3.00 


3.50 


± 


0.01±0.02 


0.03±0.02 


0.17±0.06 






3.50 


4.00 


± 


0.005±0.011 


0.020±0.014 


0.11±0.06 






4.00 


5.00 


± 




0.01±0.01 


0.04±0.04 






5.00 


6.50 


± 


± 


0.00±0.02 


0.01±0.01 






6.50 


8.00 


± 


± 


± 


0.00±0.03 
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TABLE XXIII: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 er 7r / (dpdQ,) , for p-Pb interactions at 3, 5, 8, 12 GeV/c. Each row refers to a different (p mm < p < Pmax, "min — 
< 6*max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 




d 2 a n /{dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 




(barn/(sr 


GeV/c)) 












3 GeV/c 


5 GeV/c 


8 GeV/c 


12 GeV/c 


0.050 


0.100 


0.50 


1.00 


0.06±0.06 


0.37±0.11 


0.72±0.11 


1.35±0.22 






1.00 


1.50 




0.30±0.07 


0.54±0.08 


1.37±0.18 






1.50 


2.00 


0.01±0.02 


0.15±0.05 


0.47±0.06 


0.92±0.14 






2.00 


2.50 




0.06±0.03 


0.27±0.04 


0.77±0.12 






2.50 


3.00 


± 


0.04±0.02 


0.15±0.04 


0.43±0.09 






3.00 


3.50 


± 


*_|_* 


0.05±0.02 


0.47±0.12 






3.50 


4.00 


± 


*_|_* 


0.06±0.02 


0.29±0.17 






4.00 


5.00 


± 


*_|_* 


0.034±0.012 


0.13±0.09 






5.00 


6.50 


± 


± 




0.03±0.04 






6.50 


8.00 


± 


± 


± 


0.008±0.043 


0.100 


0.150 


0.50 


1.00 


0.24±0.13 


0.43±0.11 


1.01±0.14 


2.23±0.32 






1.00 


1.50 


0.004±0.008 


0.19±0.06 


0.58±0.07 


1.39±0.18 






1.50 


2.00 


0.003±0.007 


0.05±0.02 


0.37±0.05 


0.86±0.13 






2.00 


2.50 




0.07±0.03 


0.17±0.03 


0.50±0.09 






2.50 


3.00 


± 


0.03±0.02 


0.14±0.03 


0.31±0.07 






3.00 


3.50 


± 


0.007±0.006 


0.07±0.02 


0.21±0.05 






3.50 


4.00 


± 


0.003±0.006 


0.037±0.011 


0.19±0.05 






4.00 


5.00 


± 


0.000±0.002 


0.013±0.006 


0.09±0.03 






5.00 


6.50 


± 


± 


0.002±0.002 


0.024±0.011 






6.50 


8.00 


± 


± 


± 


0.002±0.002 


0.150 


0.200 


0.50 


1.00 


0.11±0.08 


0.39±0.10 


0.98±0.13 


1.98±0.30 






1.00 


1.50 


0.07±0.06 


0.14±0.05 


0.50±0.07 


0.90±0.15 






1.50 


2.00 




0.05±0.02 


0.27±0.05 


0.58±0.11 






2.00 


2.50 




0.04±0.02 


0.17±0.03 


0.33±0.08 






2.50 


3.00 


± 


0.01±0.01 


0.07±0.02 


0.24±0.06 






3.00 


3.50 


± 


0.007±0.010 


0.033±0.012 


0.19±0.05 






3.50 


4.00 


± 


0.002±0.006 


0.011±0.007 


0.08±0.03 






4.00 


5.00 


± 






0.05±0.02 






5.00 


6.50 


± 


± 


0.000±0.001 


0.019±0.012 






6.50 


8.00 


± 


± 


± 


0.002±0.002 


0.200 


0.250 


0.50 


1.00 


0.19±0.14 


0.41±0.12 


0.68±0.11 


1.41±0.25 






1.00 


1.50 


0.03±0.04 


0.16±0.06 


0.37±0.07 


0.85±0.16 






1.50 


2.00 


0.01±0.04 


0.02±0.03 


0.17±0.05 


0.81±0.19 






2.00 


2.50 




0.02±0.02 


0.07±0.02 


0.40±0.12 






2.50 


3.00 


± 


0.021±0.026 


0.014±0.008 


0.16±0.06 






3.00 


3.50 


± 


*_|_* 


*_L_* 


0.15±0.08 






3.50 


4.00 


± 


*_|_* 


0.00±0.01 


0.04±0.04 






4.00 


5.00 


± 


*_|_=f= 


0.00±0.01 


0.01±0.02 






5.00 


6.50 


± 


± 




0.00±0.01 






6.50 


8.00 


± 


± 


± 


0.00±0.02 
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TABLE XXIV: HARP results for the double-differential ir + production cross-section in the laboratory system, 
d a w I (dpdfl) , for p-Be interactions at 8,8.9,12 GeV/c. Each row refers to a different (p mm < p < Pmax,# mm < 
< #max) bin, where p and are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 


< 


fa^ /(dpdQ) 


(rad) 


(rad) (GeV/c) (GeV/c) 


(barn/(sr GeV/c)) 










8 GeV/c 


8.9 GeV/c 


12 GeV/c 


0.025 


0.050 


0.50 


0.75 


0.13±0.04 


0.12±0.03 


0.17±0.05 






0.75 


1.00 


0.15±0.03 


0.13±0.03 


0.15±0.03 






1.00 


1.25 


0.18±0.03 


0.12±0.03 


0.14±0.03 






1.25 


1.50 


0.13±0.02 


0.13±0.03 


0.18±0.03 






1.50 


2.00 


0.20±0.02 


0.17±0.02 


0.22±0.03 






2.00 


2.50 


0.19±0.03 


0.19±0.02 


0.25±0.03 






2.50 


3.00 


0.18±0.02 


0.18±0.02 


0.25±0.03 






3.00 


3.50 


0.09±0.02 


0.15±0.02 


0.27±0.03 






3.50 


4.00 


0.047±0.014 


0.08±0.02 


0.24±0.03 






4.00 


5.00 


0.015±0.008 


0.069±0.012 


0.15±0.02 






5.00 


6.50 


0.013±0.006 


0.022±0.010 


0.083±0.010 






6.50 


8.00 


± 


± 


0.038±0.006 


0.050 


0.075 


0.50 


0.75 


0.13±0.03 


0.16±0.03 


0.16±0.03 






0.75 


1.00 


0.20±0.02 


0.14±0.02 


0.13±0.02 






1.00 


1.25 


0.14±0.02 


0.17±0.02 


0.14±0.02 






1.25 


1.50 


0.16±0.02 


0.15±0.02 


0.13±0.02 






1.50 


2.00 


0.20±0.02 


0.19±0.02 


0.21±0.02 






2.00 


2.50 


0.19±0.02 


0.22±0.02 


0.22±0.02 






2.50 


3.00 


0.15±0.02 


0.18±0.02 


0.23±0.02 






3.00 


3.50 


0.067±0.013 


0.104±0.013 


0.18±0.02 






3.50 


4.00 


0.032±0.009 


0.096±0.011 


0.16±0.02 






4.00 


5.00 


0.029±0.007 


0.050±0.008 


0.101±0.010 






5.00 


6.50 


0.004±0.001 


0.015±0.003 


0.047±0.006 






6.50 


8.00 


± 


± 


0.018±0.003 


0.075 


0.100 


0.50 


0.75 


0.15±0.02 


0.20±0.03 


0.22±0.03 






0.75 


1.00 


0.17±0.02 


0.18±0.02 


0.17±0.02 






1.00 


1.25 


0.17±0.02 


0.17±0.02 


0.19±0.02 






1.25 


1.50 


0.19±0.02 


0.17±0.02 


0.23±0.02 






1.50 


2.00 


0.21±0.02 


0.21±0.02 


0.25±0.02 






2.00 


2.50 


0.140±0.012 


0.190±0.013 


0.22±0.02 






2.50 


3.00 


0.093±0.011 


0.138±0.009 


0.196±0.015 






3.00 


3.50 


0.043±0.007 


0.077±0.008 


0.144±0.013 






3.50 


4.00 


0.025±0.004 


0.044±0.005 


0.092±0.009 






4.00 


5.00 


0.016±0.003 


0.027±0.003 


0.053±0.005 






5.00 


6.50 


0.002±0.000 


0.005±0.001 


0.023±0.003 






6.50 


8.00 


± 


± 


0.006±0.001 
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TABLE XXV: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d a 71 /(dpdSY), for p-Be interactions at 8,8.9,12 GeV/c. Each row refers to a different (p mm < p < J3max,# mm < 
< #max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



^min 


#max 


Pmin 


Pmax 


dV* /(dpdQ) 




(rad) 


(rad) (GeV/c) (GeV/c) 


(barn/(sr GeV/c)) 












8 GeV/c 


8.9 GeV/c 


12 GeV/c 


0.025 


0.050 


0.50 


0.75 


0.08±0.03 


0.11±0.02 


0.07±0 


02 






0.75 


1.00 


0.12±0.02 


0.10±0.02 


0.14±0 


03 






1.00 


1.25 


0.14±0.03 


0.17±0.03 


0.16±0 


03 






1.25 


1.50 


0.10±0.02 


0.15±0.02 


0.20±0 


04 






1.50 


2.00 


0.12±0.02 


0.11±0.02 


0.13±0 


02 






2.00 


2.50 


0.075±0.014 


0.09±0.02 


0.12±0 


02 






2.50 


3.00 


0.058±0.012 


0.074±0.013 


0.11±0 


02 






3.00 


3.50 


0.022±0.006 


0.057±0.011 


0.10±0 


02 






3.50 


4.00 


0.020±0.007 


0.042±0.008 


0.057±0 


011 






4.00 


5.00 


0.016±0.004 


0.028±0.004 


0.055±0 


008 






5.00 


6.50 


0.001±0.001 


0.006±0.002 


0.028±0 


006 






6.50 


8.00 


± 


± 


0.009±0 


003 


0.050 


0.075 


0.50 


0.75 


0.11±0.02 


0.08±0.02 


0.12±0 


03 






0.75 


1.00 


0.12±0.02 


0.11±0.02 


0.16±0 


02 






1.00 


1.25 


0.12±0.02 


0.14±0.02 


0.19±0 


03 






1.25 


1.50 


0.079±0.014 


0.126±0.015 


0.17±0 


02 






1.50 


2.00 


0.106±0.013 


0.121±0.013 


0.18±0 


02 






2.00 


2.50 


0.063±0.009 


0.092±0.009 


0.12±0 


02 






2.50 


3.00 


0.036±0.008 


0.053±0.007 


0.104±0 


014 






3.00 


3.50 


0.027±0.006 


0.045±0.006 


0.089±0 


012 






3.50 


4.00 


0.018±0.005 


0.034±0.006 


0.051±0 


008 






4.00 


5.00 


0.009±0.003 


0.021±0.003 


0.047±0 


006 






5.00 


6.50 


0.001±0.001 


0.003±0.001 


0.016±0 


003 






6.50 


8.00 


± 


± 


0.006±0 


002 


0.075 


0.100 


0.50 


0.75 


0.13±0.02 


0.11±0.02 


0.14±0 


03 






0.75 


1.00 


0.13±0.02 


0.15±0.02 


0.13±0 


02 






1.00 


1.25 


0.12±0.02 


0.162±0.015 


0.15±0 


02 






1.25 


1.50 


0.13±0.02 


0.136±0.013 


0.14±0 


02 






1.50 


2.00 


0.096±0.010 


0.116±0.009 


0.145±0 


014 






2.00 


2.50 


0.057±0.007 


0.075±0.007 


0.126±0 


013 






2.50 


3.00 


0.031±0.007 


0.046±0.005 


0.092±0 


010 






3.00 


3.50 


0.012±0.003 


0.032±0.003 


0.067±0 


008 






3.50 


4.00 


0.007±0.002 


0.018±0.004 


0.040±0 


006 






4.00 


5.00 


0.005±0.002 


0.006±0.001 


0.025±0 


003 






5.00 


6.50 






0.006±0 


001 






6.50 


8.00 


± 


± 
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TABLE XXVI: HARP results for the double-differential ir + and 7r~ production cross-section in the laboratory 
system, d 2 a n /(dpdQ), for p-C interactions at 8 and 12 GeV/c. Each row refers to a different (p m j n < p < 
Pmax,# mm < < 6<max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of the 
diagonal elements of the covariance matrix are given. 



y min y max Pmin Pmax 
(rad) (rad) (GeV/c) (GeV/c) 



d 2 a" /(dpdQ,) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 



d 2 ^' /(dpdQ) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 



0.025 


0.050 





50 





75 


0.13±0 


04 


0.22±0 


07 


0.07±0 


03 


0.14±0 


05 









75 


1 


00 


0.15±0 


04 


0.16±0 


04 


0.16±0 


04 


0.15±0 


04 






1 


00 


1 


25 


0.16±0 


04 


0.20±0 


05 


0.19±0 


04 


0.20±0 


04 






1 


25 


1 


50 


0.24±0 


04 


0.25±0 


05 


0.08±0 


02 


0.22±0 


05 






1 


50 


2 


00 


0.22±0 


03 


0.24±0 


04 


0.16±0 


03 


0.26±0 


04 






2 


00 


2 


50 


0.19±0 


03 


0.25±0 


04 


0.08±0 


02 


0.10±0 


03 






2 


50 


3 


00 


0.21±0 


03 


0.34±0 


04 


0.07±0 


02 


0.11±0 


03 






3 


00 


3 


50 


0.09±0 


02 


0.27±0 


04 


0.034±0 


010 


0.13±0 


03 






3 


50 


4 


00 


0.05±0 


02 


0.29±0 


04 


0.020±0 


009 


0.09±0 


02 






4 


00 


5 


00 


0.019±0 


011 


0.15±0 


02 


0.017±0 


005 


0.062±0 


012 






5 


00 


6 


50 


0.021±0 


010 


0.12±0 


02 


0.001±0 


001 


0.031±0 


008 






6 


50 


8 


00 


± 




0.055±0 


010 


± 




0.010±0 


005 


0.050 


0.075 





50 





75 


0.14±0 


03 


0.23±0 


05 


0.10±0 


03 


0.10±0 


03 









75 


1 


00 


0.18±0 


03 


0.21±0 


04 


0.11±0 


02 


0.16±0 


03 






1 


00 


1 


25 


0.18±0 


03 


0.18±0 


03 


0.15±0 


03 


0.21±0 


03 






1 


25 


1 


50 


0.17±0 


02 


0.22±0 


04 


0.17±0 


03 


0.20±0 


03 






1 


50 


2 


00 


0.23±0 


02 


0.25±0 


03 


0.12±0 


02 


0.17±0 


02 






2 


00 


2 


50 


0.19±0 


02 


0.29±0 


03 


0.080±0 


013 


0.15±0 


02 






2 


50 


3 


00 


0.19±0 


02 


0.28±0 


03 


0.037±0 


010 


0.12±0 


02 






3 


00 


3 


50 


0.07±0 


02 


0.20±0 


02 


0.031±0 


008 


0.11±0 


02 






3 


50 


4 


00 


0.039±0 


013 


0.19±0 


02 


0.036±0 


009 


0.075±0 


014 






4 


00 


5 


00 


0.035±0 


008 


0.105±0 


013 


0.010±0 


005 


0.053±0 


009 






5 


00 


6 


50 


0.007±0 


002 


0.055±0 


007 






0.012±0 


004 






6 


50 


8 


00 


± 




0.029±0 


005 


± 




0.004±0 


002 


0.075 


0.100 





50 





75 


0.13±0 


03 


0.15±0 


03 


0.19±0 


03 


0.14±0 


03 









75 


1 


00 


0.26±0 


03 


0.19±0 


03 


0.12±0 


02 


0.12±0 


02 






1 


00 


1 


25 


0.22±0 


02 


0.26±0 


03 


0.15±0 


02 


0.18±0 


03 






1 


25 


1 


50 


0.22±0 


02 


0.26±0 


03 


0.15±0 


02 


0.20±0 


03 






1 


50 


2 


00 


0.27±0 


02 


0.30±0 


03 


0.100±0 


013 


0.21±0 


02 






2 


00 


2 


50 


0.18±0 


02 


0.25±0 


02 


0.065±0 


010 


0.15±0 


02 






2 


50 


3 


00 


0.095±0 


013 


0.20±0 


02 


0.029±0 


008 


0.09±0 


02 






3 


00 


3 


50 


0.054±0 


010 


0.143±0 


014 


0.019±0 


004 


0.056±0 


009 






3 


50 


4 


00 


0.023±0 


005 


0.112±0 


013 


0.017±0 


004 


0.042±0 


007 
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TABLE XXVII: HARP results for the double-differential tt + production cross-section in the laboratory system, 
d 2 <r n /(dpdQ), for p-Al interactions at 8,12,12.9 GeV/c. Each row refers to a different (p mm < p < pmax,# mm < 
< 9ma,x) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



min 


#max 


Pmin 


Pmax 


d 2 a' K+ /(dpdQ) 


(rad) 


(rad) (GeV/c) 


(GeV/c) 


(barn/(sr GeV/c)) 










8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.025 


0.050 


0.50 


0.75 


0.34±0.10 


0.46±0.15 


0.28±0.08 






0.75 


1.00 


0.30±0.07 


0.44±0.10 


0.43±0.07 






1.00 


1.25 


0.23±0.06 


0.31±0.08 


0.43±0.05 






1.25 


1.50 


0.25±0.05 


0.39±0.09 


0.36±0.05 






1.50 


2.00 


0.31±0.05 


0.45±0.08 


0.45±0.04 






2.00 


2.50 


0.26±0.04 


0.50±0.09 


0.50±0.05 






2.50 


3.00 


0.30±0.05 


0.49±0.08 


0.47±0.04 






3.00 


3.50 


0.18±0.03 


0.42±0.07 


0.52±0.04 






3.50 


4.00 


0.25±0.04 


0.42±0.06 


0.42±0.04 






4.00 


5.00 


0.12±0.03 


0.27±0.04 


0.30±0.02 






5.00 


6.50 


0.027±0.013 


0.12±0.02 


0.19±0.02 






6.50 


8.00 


± 


0.06±0.02 


0.108±0.013 


0.050 


0.075 


0.50 


0.75 


0.31±0.07 


0.44±0.10 


0.48±0.08 






0.75 


1.00 


0.33±0.05 


0.33±0.06 


0.45±0.04 






1.00 


1.25 


0.24±0.04 


0.31±0.06 


0.31±0.03 






1.25 


1.50 


0.33±0.05 


0.37±0.07 


0.48±0.05 






1.50 


2.00 


0.38±0.04 


0.44±0.06 


0.47±0.04 






2.00 


2.50 


0.37±0.04 


0.46±0.06 


0.49±0.03 






2.50 


3.00 


0.27±0.03 


0.48±0.06 


0.44±0.03 






3.00 


3.50 


0.20±0.03 


0.33±0.05 


0.44±0.03 






3.50 


4.00 


0.12±0.02 


0.28±0.05 


0.34±0.02 






4.00 


5.00 


0.11±0.02 


0.18±0.03 


0.23±0.02 






5.00 


6.50 


0.009±0.003 


0.067±0.012 


0.111±0.008 






6.50 


8.00 


± 


0.035±0.008 


0.049±0.006 


0.075 


0.100 


0.50 


0.75 


0.32±0.06 


0.47±0.09 


0.55±0.06 






0.75 


1.00 


0.38±0.05 


0.57±0.08 


0.48±0.04 






1.00 


1.25 


0.40±0.04 


0.52±0.07 


0.48±0.04 






1.25 


1.50 


0.41±0.04 


0.43±0.06 


0.55±0.04 






1.50 


2.00 


0.34±0.03 


0.51±0.05 


0.56±0.03 






2.00 


2.50 


0.27±0.03 


0.48±0.05 


0.45±0.03 






2.50 


3.00 


0.20±0.02 


0.41±0.04 


0.38±0.02 






3.00 


3.50 


0.16±0.02 


0.25±0.03 


0.30±0.02 






3.50 


4.00 


0.10±0.02 


0.19±0.02 


0.23±0.02 






4.00 


5.00 


0.056±0.008 


0.117±0.015 


0.126±0.011 






5.00 


6.50 


0.004±0.002 


0.034±0.007 


0.052±0.005 






6.50 


8.00 


± 


0.010±0.003 


0.015±0.002 
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TABLE XXVIII: HARP results for the double-differential ir~ production cross-section in the laboratory system, 
d 2 a 71 / (dpdQ.) , for p-Al interactions at 8,12,12.9 GeV/c. Each row refers to a different (p mm < p < pmax,# mm < 
9 < 9max) bin, where p and 9 are the pion momentum and polar angle, respectively. The central value as well as 
the square-root of the diagonal elements of the covariance matrix are given. 



min 


#max 


^min 


Pmax 


d 2 a n /{dpdQ) 


(rad) 


(rad) (GeV/c) 


(GeV/c) 


(barn/(sr GeV/c)) 










8 GeV/c 


12 GeV/c 


12.9 GeV/c 


0.025 


0.050 


0.50 


0.75 


0.11±0.05 


0.09±0.05 


0.24±0.06 






0.75 


1.00 


0.27±0.06 


0.24±0.07 


0.36±0.05 






1.00 


1.25 


0.35±0.07 


0.43±0.10 


0.37±0.05 






1.25 


1.50 


0.20±0.05 


0.57±0.12 


0.49±0.06 






1.50 


2.00 


0.20±0.04 


0.28±0.07 


0.39±0.04 






2.00 


2.50 


0.14±0.04 


0.12±0.04 


0.27±0.04 






2.50 


3.00 


0.10±0.02 


0.15±0.05 


0.23±0.03 






3.00 


3.50 


0.09±0.02 


0.20±0.05 


0.25±0.03 






3.50 


4.00 


0.05±0.02 


0.10±0.03 


0.15±0.03 






4.00 


5.00 


0.015±0.007 


0.11±0.03 


0.103±0.014 






5.00 


6.50 




0.041±0.014 


0.065±0.009 






6.50 


8.00 


± 


0.009±0.006 


0.021±0.006 


0.050 


0.075 


0.50 


0.75 


0.19±0.05 


0.45±0.11 


0.35±0.06 






0.75 


1.00 


0.29±0.05 


0.39±0.07 


0.41±0.04 






1.00 


1.25 


0.18±0.04 


0.41±0.07 


0.42±0.04 






1.25 


1.50 


0.15±0.03 


0.22±0.05 


0.40±0.04 






1.50 


2.00 


0.20±0.03 


0.30±0.05 


0.33±0.03 






2.00 


2.50 


0.13±0.02 


0.28±0.05 


0.31±0.03 






2.50 


3.00 


0.07±0.02 


0.19±0.04 


0.22±0.02 






3.00 


3.50 


0.051±0.012 


0.16±0.03 


0.21±0.02 






3.50 


4.00 


0.047±0.011 


0.15±0.03 


0.14±0.02 






4.00 


5.00 


0.015±0.007 


0.07±0.02 


0.085±0.009 






5.00 


6.50 


0.002±0.002 


0.037±0.009 


0.043±0.005 






6.50 


8.00 


± 


0.006±0.004 


0.013±0.003 


0.075 


0.100 


0.50 


0.75 


0.32±0.06 


0.35±0.07 


0.34±0.05 






0.75 


1.00 


0.29±0.04 


0.34±0.05 


0.41±0.04 






1.00 


1.25 


0.25±0.03 


0.39±0.06 


0.38±0.03 






1.25 


1.50 


0.22±0.03 


0.35±0.05 


0.40±0.03 






1.50 


2.00 


0.18±0.02 


0.36±0.04 


0.38±0.03 






2.00 


2.50 


0.11±0.02 


0.26±0.04 


0.28±0.02 






2.50 


3.00 


0.07±0.02 


0.15±0.03 


0.18±0.02 






3.00 


3.50 


0.032±0.007 


0.13±0.02 


0.153±0.011 






3.50 


4.00 


0.021±0.005 


0.072±0.014 


0.096±0.010 






4.00 


5.00 


0.009±0.003 


0.030±0.007 


0.057±0.006 






5.00 


6.50 




0.008±0.003 


0.018±0.002 






6.50 


8.00 


± 


O.OOliO.001 


0.003±0.001 
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TABLE XXIX: HARP results for the double-differential ir + and 7r~ production cross-section in the laboratory 
system, d 2 a* / \dpdQ.) , for p-Cu interactions at 8 and 12 GeV/c. Each row refers to a different (p m i n < p < 
Pmax,# mm < < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of the 
diagonal elements of the covariance matrix are given. 



p min 
(rad) 


#max Pmin Pmax 
(rad) (GeV/c) (GeV/c) 


d 2 a n+ /(dpdQ) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 


d 2 a"~ /(dpdti) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 


0.025 


0.050 


0.50 


0.75 


0.58±0.18 


0.50±0.18 


0.23±0.10 


0.28±0.13 






0.75 


1.00 


0.40±0.10 


0.56±0.14 


0.48±0.10 


0.43±0.12 






1.00 


1.25 


0.31±0.08 


0.45±0.12 


0.43±0.10 


0.69±0.15 






1.25 


1.50 


0.47±0.10 


0.80±0.17 


0.22±0.07 


0.89±0.18 






1.50 


2.00 


0.73±0.10 


0.77±0.13 


0.35±0.07 


0.64±0.13 






2.00 


2.50 


0.34±0.07 


0.79±0.13 


0.19±0.05 


0.28±0.09 






2.50 


3.00 


0.44±0.08 


1.03±0.15 


0.20±0.04 


0.38±0.09 






3.00 


3.50 


0.18±0.05 


0.78±0.11 


0.09±0.03 


0.24±0.07 






3.50 


4.00 


0.13±0.06 


0.66±0.10 


0.05±0.03 


0.17±0.05 






4.00 


5.00 


0.04±0.03 


0.41±0.06 


0.030±0.012 


0.20±0.04 






5.00 


6.50 


0.02±0.02 


0.22±0.04 


* ± * 


0.08±0.02 






6.50 


8.00 


± 


0.12±0.03 


± 


0.016±0.010 


0.050 


0.075 


0.50 


0.75 


0.60±0.13 


0.89±0.19 


0.35±0.08 


0.51±0.14 






0.75 


1.00 


0.55±0.08 


0.70±0.12 


0.33±0.06 


0.52±0.10 






1.00 


1.25 


0.39±0.06 


0.59±0.11 


0.39±0.07 


0.70±0.12 






1.25 


1.50 


0.53±0.08 


0.74±0.13 


0.35±0.06 


0.48±0.09 






1.50 


2.00 


0.62±0.07 


0.69±0.09 


0.37±0.05 


0.67±0.10 






2.00 


2.50 


0.46±0.06 


0.64±0.09 


0.22±0.04 


0.49±0.08 






2.50 


3.00 


0.32±0.05 


0.64±0.08 


0.09±0.02 


0.35±0.06 






3.00 


3.50 


0.15±0.05 


0.53±0.07 


0.05±0.02 


0.30±0.06 






3.50 


4.00 


0.06±0.03 


0.41±0.07 


0.08±0.02 


0.18±0.04 






4.00 


5.00 


0.05±0.02 


0.29±0.04 


0.016±0.010 


0.13±0.03 






5.00 


6.50 


0.006±0.003 


0.11±0.02 


0.000±0.001 


0.044±0.014 






6.50 


8.00 


± 


0.043±0.012 


± 


0.014±0.010 


0.075 


0.100 


0.50 


0.75 


0.43±0.09 


0.65±0.12 


0.52±0.09 


0.72±0.14 






0.75 


1.00 


0.64±0.08 


0.84±0.12 


0.53±0.07 


0.53±0.08 






1.00 


1.25 


0.65±0.07 


0.83±0.10 


0.36±0.05 


0.73±0.10 






1.25 


1.50 


0.71±0.07 


0.83±0.10 


0.41±0.06 


0.82±0.10 






1.50 


2.00 


0.51±0.05 


0.90±0.09 


0.28±0.04 


0.49±0.06 






2.00 


2.50 


0.43±0.05 


0.82±0.08 


0.14±0.02 


0.35±0.05 






2.50 


3.00 


0.21±0.03 


0.59±0.06 


0.08±0.02 


0.35±0.05 






3.00 


3.50 


0.11±0.02 


0.53±0.06 


0.033±0.009 


0.17±0.03 






3.50 


4.00 


0.048±0.014 


0.32±0.04 


0.032±0.009 


0.13±0.02 



■55 



TABLE XXX: HARP results for the double-differential ir + and 7r~ production cross-section in the laboratory 
system, d 2 a w / (dpdQ.) , for p- Sn interactions at 8 and 12 GeV/c. Each row refers to a different (p mm < p < 
Pmax,# mm < < fciax) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of the 
diagonal elements of the covariance matrix are given. 



^min fc'max p m i n Pmax 
(rad) (rad) (GeV/c) (GeV/c) 



0.025 0.050 



0.050 0.075 



0.075 0.100 



dV* /(dpdQ) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 



d 2 cr n ~ /(dpdQ) 
(barn/(sr GeV/c)) 
8 GeV/c 12 GeV/c 



0.50 





75 


1.06±0 


27 


1 


06±0 


25 


0.86±0.24 


1.03±0 


24 


0.75 


1 


00 


0.89±0 


18 





93±0 


17 


0.81±0.15 


0.79±0 


15 


1.00 


1 


25 


0.73±0 


16 





72±0 


14 


0.69±0.15 


1.22±0 


19 


1.25 


1 


50 


0.51±0 


13 


1 


10±0 


18 


0.37±0.11 


1.15±0 


18 


1.50 


2 


00 


0.62±0 


11 


1 


00±0 


13 


0.51±0.10 


0.98±0 


15 


2.00 


2 


50 


0.56±0 


11 


1 


00±0 


13 


0.35±0.08 


0.33±0 


09 


2.50 


3 


00 


0.47±0 


10 


1 


15±0 


15 


0.24±0.06 


0.46±0 


09 


3.00 


3 


50 


0.30±0 


09 





87±0 


11 


0.15±0.04 


0.40±0 


09 


3.50 


4 


00 


0.23±0 


27 





73±0 


10 


0.05±0.04 


0.31±0 


06 


4.00 


5 


00 


0.06±0 


09 





64±0 


07 


0.04±0.02 


0.23±0 


04 


5.00 


6 


50 


0.07±0 


05 





27±0 


05 


0.00±0.01 


0.09±0 


02 


6.50 


8 


00 


± 







14±0 


03 


± 


0.022±0 


012 


0.50 





75 


0.82±0 


17 


1 


58±0 


24 


0.51±0.12 


0.79±0 


15 


0.75 


1 


00 


0.76±0 


11 





93±0 


13 


0.64±0.11 


0.81±0 


11 


1.00 


1 


25 


0.82±0 


11 





72±0 


11 


0.49±0.10 


1.17±0 


14 


1.25 


1 


50 


0.84±0 


12 





99±0 


14 


0.38±0.08 


0.69±0 


10 


1.50 


2 


00 


0.73±0 


09 





91±0 


09 


0.49±0.07 


0.84±0 


10 


2.00 


2 


50 


0.70±0 


08 





99±0 


10 


0.34±0.05 


0.56±0 


08 


2.50 


3 


00 


0.52±0 


08 





95±0 


09 


0.16±0.04 


0.46±0 


06 


3.00 


3 


50 


0.21±0 


07 





94±0 


09 


0.13±0.03 


0.41±0 


06 


3.50 


4 


00 


0.05±0 


04 





53±0 


08 


0.06±0.02 


0.29±0 


04 


4.00 


5 


00 


0.09±0 


03 





40±0 


04 


0.03±0.02 


0.17±0 


03 


5.00 


6 


50 


0.011±0 


005 





18±0 


03 


0.006±0.005 


0.044±0 


012 


6.50 


8 


00 


± 







08±0 


02 


± 


0.014±0 


006 


0.50 





75 


0.72±0 


13 


1 


05±0 


15 


0.49±0.10 


0.94±0 


15 


0.75 


1 


00 


0.89±0 


11 


1 


21±0 


13 


0.75±0.09 


1.06±0 


12 


1.00 


1 


25 


0.72±0 


09 


1 


25±0 


12 


0.53±0.08 


1.15±0 


12 


1.25 


1 


50 


0.95±0 


11 


1 


19±0 


12 


0.56±0.08 


1.03±0 


10 


1.50 


2 


00 


0.75±0 


07 


1 


21±0 


09 


0.36±0.05 


0.80±0 


07 


2.00 


2 


50 


0.66±0 


07 


1 


04±0 


08 


0.17±0.03 


0.46±0 


05 


2.50 


3 


00 


0.28±0 


06 





84±0 


07 


0.11±0.02 


0.46±0 


05 


3.00 


3 


50 


0.09±0 


02 





54±0 


05 


0.06±0.02 


0.22±0 


03 


3.50 


4 


00 


0.06±0 


02 





38±0 


04 


0.028±0.009 


0.14±0 


02 
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TABLE XXXI: HARP results for the double-differential ir + and 7r~ production cross-section in the laboratory 
system, d 2 a™ / (dpdQ.) , for p-Ta interactions at 8 and 12 GeV/c. Each row refers to a different (p mm < p < 
Pmax^jjjin < < #max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of the 
diagonal elements of the covariance matrix are given. 



n 

a min 


#max 


Pmin 


Pmax 


d 2 a w+ /(dpdQ) 


d 2 a* 


/(dpdQ) 


(rad) 


(rad) (GeV/c) (GeV/c) 


(barn/(sr 


GeV/c)) 


(barn/(sr GeV/c)) 










8 GeV/c 


12 GeV/c 


8 GeV/c 


12 GeV/c 


0.025 


0.050 


0.50 


0.75 


1.39±0.40 


1.74±0.48 


0.69±0.27 


0.94±0.30 






0.75 


1.00 


0.71±0.20 


1.44±0.32 


0.65±0.17 


1.19±0.28 






1.00 


1.25 


0.57±0.18 


1.07±0.26 


0.52±0.16 


1.98±0.37 






1.25 


1.50 


0.60±0.17 


1.28±0.28 


0.44±0.15 


1.85±0.36 






1.50 


2.00 


0.92±0.17 


1.27±0.23 


0.61±0.13 


0.73±0.19 






2.00 


2.50 


0.60±0.14 


1.18±0.22 


0.28±0.09 


0.53±0.17 






2.50 


3.00 


0.63±0.15 


1.47±0.24 


0.22±0.07 


0.47±0.14 






3.00 


3.50 


0.25±0.13 


1.09±0.19 


0.17±0.06 


0.48±0.14 






3.50 


4.00 


0.02±0.05 


1.04±0.18 


0.05±0.05 


0.35±0.10 






4.00 


5.00 


0.01±0.24 


0.71±0.12 


0.03±0.04 


0.27±0.07 






5.00 


6.50 


0.02±0.07 


0.32±0.07 


0.00±0.03 


0.11±0.04 






6.50 


8.00 


± 


0.17±0.05 


± 


0.02±0.02 


0.050 


0.075 


0.50 


0.75 


0.82±0.22 


1.89±0.37 


0.61±0.17 


0.85±0.22 






0.75 


1.00 


0.99±0.17 


1.26±0.23 


0.77±0.15 


1.61±0.26 






1.00 


1.25 


0.76±0.13 


0.88±0.17 


0.79±0.16 


1.87±0.26 






1.25 


1.50 


0.90±0.15 


1.14±0.23 


0.51±0.12 


0.99±0.18 






1.50 


2.00 


0.84±0.12 


1.14±0.16 


0.48±0.09 


0.95±0.15 






2.00 


2.50 


0.71±0.11 


1.05±0.15 


0.35±0.07 


0.91±0.15 






2.50 


3.00 


0.44±0.09 


1.27±0.18 


0.13±0.05 


0.59±0.12 






3.00 


3.50 


0.18±0.09 


0.90±0.13 


0.05±0.02 


0.49±0.10 






3.50 


4.00 


0.06±0.06 


0.72±0.12 


0.13±0.04 


0.28±0.09 






4.00 


5.00 


0.07±0.05 


0.36±0.07 


0.01±0.02 


0.14±0.04 






5.00 


6.50 


0.004±0.004 


0.20±0.04 


*_|_* 


0.07±0.02 






6.50 


8.00 


± 


0.07±0.02 


± 


0.021±0.013 


0.075 


0.100 


0.50 


0.75 


0.55±0.14 


2.01±0.30 


0.70±0.15 


1.49±0.26 






0.75 


1.00 


1.12±0.16 


1.41±0.20 


0.64±0.10 


0.89±0.14 






1.00 


1.25 


0.96±0.13 


1.83±0.22 


0.72±0.12 


1.16±0.17 






1.25 


1.50 


1.10±0.14 


1.07±0.15 


0.41±0.08 


1.01±0.15 






1.50 


2.00 


0.77±0.09 


1.53±0.15 


0.38±0.06 


0.84±0.11 






2.00 


2.50 


0.49±0.08 


1.06±0.13 


0.21±0.04 


0.93±0.12 






2.50 


3.00 


0.27±0.06 


1.13±0.12 


0.11±0.03 


0.39±0.08 






3.00 


3.50 


0.13±0.05 


0.50±0.08 


0.04±0.02 


0.38±0.06 






3.50 


4.00 


0.05±0.02 


0.45±0.07 


0.06±0.02 


0.21±0.04 
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TABLE XXXII: HARP results for the double-differential 7r + and 7r~ production cross-section in the laboratory 
system, d 2 a 7 ' / (dpdQ) , for p-Pb interactions at 8 and 12 GeV/c. Each row refers to a different (p m j n < p < 
Pma.x,8 m [ n < < 6<max) bin, where p and 9 are the pion momentum and polar angle, respectively. A finer 
angular binning than in the previous set of tables is used. The central value as well as the square-root of the 
diagonal elements of the covariance matrix are given. 



n 


#max 


Pmin 


Pmax 


d 2 a 7T+ /(dpdQ) 


d 2 a* /(dpdQ,) 


(rad) 


(rad) (GeV/c) (GeV/c) 


(barn/(sr 


GeV/c)) 


(barn/(sr 


GeV/c)) 










8 GeV/c 


12 GeV/c 


8 GeV/c 


12 GeV/c 


0.025 


0.050 


0.50 


0.75 


1.12±0.36 


1.66±0.56 


0.23±0.15 


1.47±0.51 






0.75 


1.00 


1.17±0.28 


2.34±0.57 


1.23±0.27 


2.20±0.54 






1.00 


1.25 


0.38±0.14 


1.31±0.39 


0.81±0.21 


1.00±0.32 






1.25 


1.50 


0.62±0.19 


1.23±0.37 


0.35±0.13 


2.00±0.54 






1.50 


2.00 


1.13±0.20 


1.59±0.36 


0.54±0.13 


1.18±0.33 






2.00 


2.50 


0.54±0.15 


1.17±0.30 


0.34±0.12 


0.55±0.23 






2.50 


3.00 


0.50±0.14 


1.39±0.34 


0.16±0.06 


0.64±0.23 






3.00 


3.50 


0.33±0.13 


1.69±0.35 


0.16±0.07 


0.61±0.24 






3.50 


4.00 


0.08±0.09 


0.73±0.20 


0.03±0.04 


0.69±0.22 






4.00 


5.00 


0.05±0.19 


0.66±0.15 


0.07±0.03 


0.44±0.14 






5.00 


6.50 


*_j_* 


0.41±0.11 


0.00±0.01 


0.16±0.22 






6.50 


8.00 


± 


0.14±0.06 


± 


0.02±0.15 


0.050 


0.075 


0.50 


0.75 


0.88±0.23 


2.65±0.60 


0.47±0.15 


1.08±0.33 






0.75 


1.00 


1.14±0.18 


1.61±0.34 


0.73±0.15 


1.25±0.28 






1.00 


1.25 


0.90±0.15 


1.06±0.26 


0.58±0.14 


1.43±0.31 






1.25 


1.50 


0.88±0.16 


1.33±0.32 


0.56±0.13 


1.19±0.29 






1.50 


2.00 


0.81±0.12 


1.26±0.23 


0.60±0.10 


0.86±0.19 






2.00 


2.50 


0.67±0.11 


1.30±0.23 


0.33±0.07 


0.58±0.16 






2.50 


3.00 


0.60±0.11 


1.11±0.22 


0.18±0.06 


0.50±0.16 






3.00 


3.50 


0.28±0.11 


0.87±0.18 


0.07±0.03 


0.44±0.13 






3.50 


4.00 


0.03±0.13 


0.61±0.14 


0.10±0.04 


0.25±0.20 






4.00 


5.00 


0.12±0.06 


0.52±0.11 


0.02±0.02 


0.19±0.16 






5.00 


6.50 


0.01±0.03 


0.20±0.06 




0.07±0.09 






6.50 


8.00 


± 


0.03±0.02 


± 


0.02±0.10 


0.075 


0.100 


0.50 


0.75 


0.72±0.17 


1.50±0.32 


0.93±0.18 


1.53±0.35 






0.75 


1.00 


1.21±0.16 


1.57±0.29 


0.70±0.11 


1.44±0.26 






1.00 


1.25 


0.75±0.11 


1.51±0.27 


0.41±0.09 


1.60±0.28 






1.25 


1.50 


0.92±0.13 


0.93±0.19 


0.64±0.11 


1.21±0.22 






1.50 


2.00 


0.81±0.10 


1.62±0.24 


0.37±0.07 


0.96±0.16 






2.00 


2.50 


0.66±0.10 


1.48±0.20 


0.23±0.05 


0.91±0.16 






2.50 


3.00 


0.28±0.07 


0.96±0.14 


0.12±0.04 


0.38±0.09 






3.00 


3.50 


0.06±0.03 


0.80±0.13 


0.037±0.015 


0.50±0.19 






3.50 


4.00 


0.027±0.015 


0.52±0.10 


0.028±0.011 


0.31±0.16 
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